Statistics of Avalanches with Relaxation, and Barkhausen Noise: A Solvable Model 



m 

o 



CZ5 



cd 



I 

o 
o 



> 

On 

(N 
l> 

^' 

O 
en 

> 



Alexander DobrinevskiQ Pierre Le Doussal, and Kay Jorg Wiese 
CNRS-Laboratoire de Physique Theorique de I'Ecole Normale Superieure, 24 rue Lhomond, 75005 Paris, France 

(Dated: April 29, 2013) 

We study a generalization of the Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model of a par- 
ticle in a Brownian force landscape, including retardation effects. We show that under monotonous 
driving the particle moves forward at all times, as it does in absence of retardation (Middleton's 
theorem). This remarkable property allows us to develop an analytical treatment. The model 
with an exponentially decaying memory kernel is realized in Barkhausen experiments with eddy- 
current relaxation, and has previously been shown numerically to account for the experimentally 
observed asymmetry of Barkhausen-pulse shapes. We elucidate another qualitatively new feature: 
the breakup of each avalanche of the standard ABBM model into a cluster of sub-avalanches, sharply 
delimited for slow relaxation under quasi-static driving. These conditions are typical for earthquake 
dynamics. With relaxation and aftershock clustering, the present model includes important ingre- 
dients for an effective description of earthquakes. We analyze quantitatively the limits of slow and 
fast relaxation for stationary driving with velocity t; > 0. The n-dependent power-law exponent 
for small velocities, and the critical driving velocity at which the particle velocity never vanishes, 
are modified. We also analyze non-stationary avalanches following a step in the driving magnetic 
field. Analytically, we obtain the mean avalanche shape at fixed size, the duration distribution of 
the first sub-avalanche, and the time dependence of the mean velocity. We propose to study these 
observables in experiments, allowing to directly measure the shape of the memory kernel, and to 
trace eddy current relaxation in Barkhausen noise. 

PACS numbers: 02.50.Ey, 05.40. Jc, 75.60.Ej 



I. INTRODUCTION AND MODEL 
A. Barkhausen noise 

The Barkhausen noise p^ is a characteristic magnetic 
signal emitted when a soft magnet is slowly magne- 
tized. It can be measured and made audible as crack- 
ling through an induction coil; periods of quiescence fol- 
lowed by pulses, or avalanches, of random strength and 
duration. The statistics of the emitted signal depends 
on material properties and its state. By analyzing the 
Barkhausen signal, one can deduce for example residual 
stresses [21 [3] or grain sizes [U [5] in metallic materials. 
Understanding how particular details of the Barkhausen 
noise statistics depend on microscopic material proper- 
ties is important for such applications. 

On the other hand, Barkhausen noise pulses are just 
one example for avalanches in disordered media. Such 
avalanches also occur in the propagation of cracks dur- 
ing fracture [BHS]; in the motion of fluid contact lines 
on a rough surface [9lfl2]. and as earthquakes driven by 
motion of tectonic plates [T3HISI- Some features of the 
avalanche statistics, like size and duration distributions 
[TTI lis] , are universal for many of these phenomena [19] . 
Barkhausen noise is easily measurable experimentally, 
and provides a good way to study aspects of avalanche 
dynamics common to all these systems. 

A first advance in the theoretical description of 
Barkhausen noise was the stochastic model postulated 
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by Alessandro, Beatrice, Bertotti and Montorsi [201 121] 
{ABBM model). They proposed modeling the domain- 
wall position u{t) through the stochastic differential 
equation (SDE) 



ru{t) = 2Is [H{t) - ku{t) + F{u{t))\ . 



(1) 



We follow here the conventions of [22] and [23]. Is is 
the saturation magnetization, and H{t) the external field 
which drives the domain-wall motion. A typical choice is 
a constant ramp rate c, H{t) ~ ct ^ kvt, which leads to 
a constant average domain- wall velocity v = c/k [20\ k 
is the demagnetizing factor characterizing the strength of 
the demagnetizing field —ku generated by effective free 
magnetic charges on the sample boundary f20l [21] • The 
domain- wall motion induces a voltage proportional to its 
velocity u(t), which is the measured Barkhausen noise 
signal. Here F{u{t)) is a random local pinning force. It 
is assumed to be a Brownian motion, i.e. Gaussian with 
correlations 



[F{u) - F{u')] 



2a\u- 



This choice may seem unnatural, since the physical dis- 
order does not exhibit such long-range correlations. It 
is only recently that it has been shown [T71 [THl 121] that 
the "ABBM guess" emerges as an effective disorder to 
describe the avalanche motion of the center-of-mass of 
the interface, denoted u(t), in the mean-field limit of 
the field theory of an elastic interface with d internal di- 
mensions. This correspondence holds both for interfaces 
driven quasi-statically [18, 25 , and for static interfaces 
at zero temperature [17]. The mean-field description is 
accurate above a certain critical internal dimension dc. 



For d < dc, a systematic expansion in e = dc — d using 
the functional renormalization group yields universal cor- 
rections to the scaling exponents pS^BB] and avalanche 
size [HI [IH] and duration [151 [13 [22] distributions. 

For the particular case of magnetic domain walls, the 
predictions of the ABBM model are well verified experi- 
mentally in certain ferromagnetic materials, for example 
FeSi alloys [lU [301 [31] • These are characterized by long- 
range dipolar forces decaying as 1/r^ between parts of 
the domain wall a distance r apart. This leads [S^ to 
a critical dimension dc = 2 coinciding with the physical 
dimension of the domain wall. In this kind of systems, 
as expected, the mean-field approximation is reasonably 
well satisfied. Measurements on other types of ferromag- 
nets, for example FeCoB alloys [3T] indicate a univer- 
sality class different from the mean-field ABBM model. 
This may be explained by short-range elasticity, and a 
critical dimension dc > 2. To describe even the center 
of mass mode in this class of domain walls, one needs 
to take into account the spatial structure of the domain 
wall. Predictions for roughness exponents [271 [21] and 
avalanche statistics \XljlSj2^22^ for this non-mean- field 
universality class have been obtained using the functional 
renormalization group. 

On the other hand, even for magnets in the mean-field 
universality class, a careful measurement of Barkhausen 
pulse shapes [ini [2H [331 [31] shows that they differ from 
the simple symmetric shape predicted by the ABBM 
model [251 I3S] ■ This hints at a more complicated equa- 
tion of motion than the first-order overdamped dynamics 
usually considered for elastic interfaces in disordered me- 
dia. 

In a physical interface, there may be additional degrees 
of freedom. One example was studied in [36l[37]. Other 
examples include deformations of a plastic medium, or 
eddy currents arising during the motion of a magnetic 
domain wall. For viscoelastic media, these can be mod- 
eled by a memory term which is non-local in time [38[l39j. 
At mean-field level, this is equivalent to a model with dy- 
namical stress overshoots [40] . Such memory terms may 
lead to interesting new phenomena, like coexistence of 
pinned and moving states 38. 39, 41j . A similar mem- 
ory term, non-local in time, is argued in [35] to describe 
the dissipation of eddy currents in magnetic domain- wall 
dynamics. 



^ ' dsj{t-s)u{s)^2I,[H{t)-ku{t)+F{u{t))]. 

(2) 

The response function f, derived by solving the Maxwell 
equations in a rectangular sample [22HM1 W2\ , is 
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;i„ = (2m + \fuj,, + (2n + \fu:^, 
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They depend on the sample width a, thickness h, perme- 
ability /i and conductivity a. (|2| and ([3]) correspond to 
Eqs. (13), (17) and (21) in [23]; we refer the reader there 
for details of the derivation. 

Zapperi et al. [22^ showed numerically that avalanche 
shapes in the model ^ are asymmetric. They concluded 
that eddy-current relaxation may be one way of explain- 
ing the experimentally observed skewness of Barkhausen 
noise pulses. They also argue that similar relaxation ef- 
fects may be relevant for other physical situations where 
asymmetric pulse shapes are observecj^ 

A simplification of Eq. ([3|) occurs when considering 
only the leading contributions for small and large relax- 
ation timea^ Then one obtains [22] a natural generaliza- 
tion of the^BBM equation Q: 



Tu{t) 



ro 
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dse^(*-'')/^u(s) 



= 2/s [H{t) ~ ku{t) + F{u{t))\ . 



(4) 



Here, r is the longest relaxation time of the eddy-current 

_ ^iu f 1 j_ 1 \-i 
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modes, r 



ing coefficients given in [22] 



r and Fq are damp- 



B. The ABBM model with retardation 

For the remainder of this work, we adopt the conven- 
tions used in the study of elastic interfaces. Let us intro- 
duce a more general model than Q, 

r]u{t)+a / dsf{t-s)u{s) = F{u{t))+m'^[w{t)~u{t)]. 

J — OO 

(5) 
which describes a particle driven in a force landscape 
F{u), with retardation. At this stage F{u) is arbitrary. 
Here f{t) is a general memory kernel with the following 
properties: 

1. /(O) = 1 (without loss of generality, since a con- 
stant may be absorbed into the parameter a). 



^ For example, 1221 mentions slip velocity profiles during earth- 
quakes [141 I43j . However, it is not clear if there are physical 
reasons to expect a relaxation of the form |2l . 



We 



have /(O) <x E 



m.,n (2n + l)' 



and Jq°° f{t)dt <x 
const. Thus, for small times, 



f{t) is well approximated by ^5{t) for some constant F. On the 
other hand, for long times, only the mode that relaxes slowest 



are relaxation times for the individual eddy current 



remains. Hence, for long times one can set /(t) ^ 



-'/■r-o.o 



■^0,0 



2. f{x) — 7^ as a; — >• oo. 

3- f'{x) < for all x, i.e. memory of the past trajec- 
tory always decays with time. 

This model possesses a remarkable property for any such 
kernel f{t) and any landscape F{u). It has monotonicity, 
i.e. it satisfies the Middleton theorem: For non-negative 
driving w > 0, after an initial transient period one has 
M > at all times. A more precise statement and a 
proof are given in Appendix [A] It has very important 
consequences, both in the driven regime, and in the limit 
of quasi-static driving, i.e. small w — >■ 0^. In that limit 
it converges to the quasi-static process u{t) — > u(w(i)), 
where u{'w) is the (forward) Middleton metastable state, 
defined as the smallest (leftmost) root of 



771 u 



Fiu) 



m w 



<^ 



u{w) 



(6) 



It is independent of the precise form of the kernel f{t). 
Hence the domain- wall position u{t) is uniquely deter- 
mined by the value of the driving field w(t), due to the 
monotonicity property |44) . This process u{'w) exhibits 
jumps at a set of individual points, the avalanche loca- 
tions Wi, and the quasi-static avalanche sizes 



uiwt) " "K ) 



(7) 



are thus independent of the retardation kernel. What 
depends on the kernel is the dynamics within these 
avalanches, and that is studied here. The quasi-static 
avalanche sizes Si have a well-defined distribution P{S) 
which has been computed for a particle in various force 
landscapes [JSJ US] and for the non-trivial case of a d- 
dimensional elastic interface using functional RG meth- 
ods 29, 47, 48J. As long as the dynamics obeys the Mid- 
dleton theorem, the avalanche- size distribution remains 
independent of the details of the dynamics |18j . 

While monotonicity holds for any F(u), in this arti- 
cle we focus on the case of the Brownian force landscape 
which can be solved analytically. As in the standard 
ABBM model, we choose the effective random pinning 
force F{u) to be a random walk, i.e. Gaussian with corre- 
lator given by ([2|. [^ We call this the ABBM model with 
retardation. In view of the application to Barkhausen 
noise, the parameter a > describes the overall strength 
of the force exerted by eddy currents on the domain wall. 
For a = 0, ([5]) reduces to the equation of the standard 
ABBM model in the conventions of [JSJ US] . 

The retarded ABBM model is particularly interesting 
in view of the monotonicity property. Other ways of gen- 
eralizing the ABBM model to include inertia, e.g. by a 
second-order derivative [351; do not inherit this property 



^ It can be realized as a stationary landscape, F{u)F{u') = Aq — 
a\u — u'\ with a cutoff at scale u ~ Aq/o", or by F{u)F{u') = 
2(Tniin(u, u') (non-stationary landscape with -F(O) = 0). In both 
cases F'(u) is a white noise, and that is the important feature. 



from the standard ABBM model. This makes the ABBM 
model with retardation very special, and it will be im- 
portant for its solution in section [111] 

When considering the particularly interesting case of 
exponential relaxation motivated in [22] , we set 



/w 



-t/r 



(8) 



T is the longest time scale of eddy-current relaxation, as 
discussed above. In this approximation, ([5]) can be re- 
written as two coupled, local equations for the domain- 
wall velocity u{t), and the eddy-current pressure /i(t), 

1 ft 

h{t) = - / dse-(*-'*)/^u(s) (9) 

T J-00 

r]u{t) + aTh{t) = F{u{t)) + m^ [w{t) - u{t)] (10) 
Tdth{t)^u{t)-h{t). (11) 

Although most of our quantitative results will be derived 
for this special case only, most qualitative features carry 
over to more general kernels with sufficiently fast decay. 
By rescaling u, w and t in Eq. ((sl) (for details, see 
section III A I, one finds the characteristic time scale 



Tm = rj/m'^ and length scale Sm = o jnv' of the stan- 
dard ABBM model (a — 0). They set the scales for 
the durations and sizes of the largest avalanches. There 
are of course avalanches of smaller size (up to some mi- 
croscopic cutoff if one defines it). The velocity scale is 
Vm = o l(r\w?^ and one can define a renewal time for the 
large avalanches as t^ = Sm/v, the limit of quasi-static 
driving being r„j ^ r^,, equivalent to v/vm <C 1. In the 
retarded ABBM model ((sl) one introduces an additional 
memory time scale r and various regimes will emerge de- 
pending on how T compares with the other time scales 
(whose meaning will be changed). 

Eq. ( |11[ ) then describes a depinning model with relax- 
ation, i.e. one can think of the disorder landscape as re- 
laxing via the additional degree of freedom h{t). This is 
a feature of interest for earthquake models as discussed 
below. In this context one considers the limit of well 
separated time scales, r^ <C r <C r^. 

Other features of Barkhausen noise predicted for the 
ABBM model with retardation are quite different from 
those of the standard ABBM model. Zapperi et al. 
P^ already realized that the inclusion of eddy currents 
leads to a skewness in the avalanche shape. In this arti- 
cle, we go further and discuss changes in the avalanche 
statistics. The relaxation of eddy currents introduces an 
additional slow time scale into the model. This leads 
to avalanches which stretch further in time. In par- 
ticular, avalanches following a kick (or, more generally, 
stopped driving) never terminate, by contrast with the 
standard ABBM model. This is because of the exponen- 
tially decaying retardation kernel, which never vanishesrl 



* For a model such that /(t) = for t > to, avalanches would 
remain of finite duration. 



Avalanche sizes however, are not changed by retardation 
in the limit of quasi-static driving, as discussed above. In 
that limit, retardation leads to a break-up of avalanches 
into sub- avalanches, which can also be called aftershocks. 
Avalanches at continuous driving overlap stronger, and 
the velocity threshold for the infinite avalanche (i.e. the 
velocity u no longer vanishes) is decreased. We now de- 
scribe these effects in detail and formulate more precise 
statements. 



C. Protocols 

Let us first review qualitatively the main situations 
that we will study, and define the terminology. 

(i) stationary driving: The driving velocity is con- 
stant, w{t) = vt, and the distributions of the domain- 
wall velocity ii and of the eddy-current pressure h reach 
a steady state, which we study. If v is large enough the 
velocity will never vanish and one has a single infinite 
avalanche, also called "continuous motion". At smaller 
f > the velocity will sometimes vanish. That de- 
fines steady state avalanches. These are more properly 
called sub-avalanches of the infinite avalanche since at fi- 
nite w > they immediately restart. Only in the limit 
f = 0^ they become well separated in time and can then 
be called steady state avalanches. 

(ii) Avalanches following a kick: We consider an initial 
condition dX t — {) prepared to lie in the "Middleton 
attractor" at u = u{w{t < 0)), as discussed above. It 
can be obtained by driving the system monotonously in 
the far past with w > 0, until memory of the initial 
condition is erased; then let it relax for a long time with 
w = until time t = 0. Hence the initial condition 
is u{t = 0) = h{t = 0) = 0. At i = 0, one changes 
the external magnetic field instantaneously by wq, i.e. 
sets w{t) = ■WQ6{t). For t > 0, the external field does 
not change anymore, thus a kick in the driving velocity 
corresponds to a step in the applied force. At i = oo the 
system has settled again into the Middleton attractor 
at u(t — oo) — u{w -t- Wo), because of the properties 
discussed above. One can thus consider the total motion 
to define a single avalanche following a kick, which is thus 
unambiguously defined. The total size S = L u{t) dt 
is the same as in the absence of retardation. We will 
ask about the total duration (which becomes infinite) 
and whether the velocity has vanished at intermediate 
times, i.e. whether the avalanche has broken into sub- 
avalanches. 

Avalanches following a kick are called non- stationary 
avalanches (since driving is non-stationary) . However, in 
the limit of wq — > 0"'' they become identical to the steady- 
state avalanches obtained by stationary driving discussed 
above (conditioned to start at i = 0). 



D. Organization of this article 

The remainder of this article is structured as follows: 
In section [Hj we discuss in more detail the phenomenol- 
ogy and the qualitative physics of the ABBM model with 
retardation. We discuss the splitting of a quasi-static 
avalanche into sub-avalanches, and the effects of retar- 
dation on the stationary and the non-stationary dynam- 
ics. 



In section HI we explain how the probability distribu- 
tion of observables linear in the domain- wall velocity can 
be computed by solving a non-linear, non-local "instan- 
ton" equation. By this, the stochastic model is mapped 
onto a purely deterministic problem of non-linear dynam- 
ics. This is a generalization of the method developed in 
PSI H5] for the standard ABBM model with arbitrary 
driving. 

Section [TV] discusses how the explicit form of the mem- 
ory kernel f{t) can be extracted in an experiment from 
the response to a kick. 

Section Iv] is devoted to an analysis of the instanton 
equations in the limit ^ <c r. This means that eddy 
currents relax much more slowly than the domain wall 
moves. In this limit, we obtain the stationary distribu- 
tions of the eddy-current pressure and domain- wall veloc- 
ity, as well as their behavior following an instantaneous 
kick in the driving field. The instanton solution reflects 
the two time scales in the problem: A short time scale, 
on which eddy currents build up but do not affect the 
dynamics, and a long time scale, on which they relax 
quasi-statically. We prove that, even after the driving has 
stopped, the velocity never becomes zero permanently. 
^ we discuss the fast-relaxation limit —■ 



VI 



» 



In section 

T. In this limit, eddy currents relax much faster than 
the domain wall moves. The instanton solutions again 
exhibit two time scales, but now eddy currents are irrel- 
evant for the long-time asymptotics. Qualitative results 
(like the fact that the domain-wall motion never stops en- 
tirely) are in agreement with those for the slow-relaxation 
limit, considered in section \V\ 

In section [VlI| we discuss non-stationary avalanches fol- 
lowing an instantaneous kick in the driving. In particu- 
lar, we compute their average shape at fixed size. 

In section [VHH we show how to include an absorbing 
boundary in the instanton solution of section [Hi] This is 
required for treating avalanches during stationary driv- 
ing. We then derive the distribution of avalanche du- 
rations in the standard ABBM model at finite driving 
velocity, w > 0, and the leading corrections for weak re- 
laxation and T = Tm- We also show numerical results 
for more general situations, and give some conjectures 
on the modification of size and duration exponents by 
retardation effects. 

Last, in section |IX[ we summarize our results. We 
discuss how they can be used to learn more about the 
dependence of Barkhausen noise on eddy current dissi- 
pation. 




(c) ABBM model with retardation, r = 3. 




(d) ABBM model with retardation, r 
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Figure 1. Splitting of an avalanche into sub-avalanches through the retardation mechanism. We have set rri^ = 1 and a = 1 
and we vary the relaxation time r. 



II. PHYSICS OF THE MODEL AND SUMMARY 
OF THE RESULTS 



Quasi-static driving: Sub-avalanches and 
aftershocks 



Consider the system either under stationary driving at 
w = 0"*", or following an infinitesimal kick wq — 0"*" as 
discussed above, and call t = the starting time of the 
avalanche. The main physics can be understood from 
figure [l] and keeping in mind the equations (11). 

In Fig.fTk we represent the usual construction for u(w) 
in the standard ABBM as the left-most solution of the 
equation ([6]) (in the figure we set m = 1). Assuming 
T,n <C Ty this construction indicates the position of the 
domain wall as a function oi w = w(t) on time scales of 
order r„. At wi the solution jumps from ui = w(wj~) to 
u'l — u{w^) ^ ui + S corresponding to an avalanche 
of size S; the latter occurs on the much faster time 
scale Tm- During the avalanche the velocity u{u) (set- 
ting 77 = 1) is given by the difference in height between 
the line m'^w = m^wi and the landscape rn^u — F{u), 
providing a graphic representation of the motion. The 
velocity iiiu) vanishes ai u = ui and u = u'^. For illus- 
tration we have represented a force landscape which is 



ABBM like at large scales but smooth at small scales. 
For the continuous ABBM model the construction is re- 
peated at all scales and one has avalanches of all smaller 
sizes. 

Let us now add retardation, setting a > 0, and varying 
the memory time r. The graphical construction corre- 
sponding to Eq. (fTTl) is represented in Fig.IIJD toflH. The 
difference in height is now the sum of ii and arh (in the 
Figure we chose a = 1), which evolves according to the 
second equation in Eq. (11). It can be rewritten as 



rduh = 1 



u{u) 



(12) 



Hence h increases from h — Q^ initially as /i ~ (u — ui)/t 
(since ii ~ ^u — ui). Thus the curve w — arh versus u 
starts with a negative slope —a. 

Another way to see this is to note that for i ^ r, the 
second equation of (11) gives 



Th{t) 



u{t) + 0{t/T) =u(t) 



Oit/r). (13) 



Inserting this into the first equation of ( |11[ ), we obtain 

riu{t) = F{u{t)) + m^ [w{t) - u{t)] - a [u{t) - Ui] + .... 

(14) 



Effectively, for short times the mass is modified from 
m^ — >■ m^ + a. Thus, while w is fixed, the end of 
the first sub-avalanche is determined not by the roots 
of rri^w = m?u — F(v), but by the roots of m?w — 
{vn? + a)u — aui — F(u). Equivalently, in the landscape 
rri^u — F{u), instead of looking at intersections with the 
horizontal curve m?w, we should look at intersections 
with m^w — a{u — Mi), a line with slope —a. 

At the point where this curve intersects first the land- 
scape m'^u — F{u) we get a point Usi < u'l where u first 
vanishes. This defines the size Si = Usi — wi of the first 
sub-avalanche. If r is small this usually occurs near the 
end, but if r is larger the original avalanche (called main 
avalanche) is divided - in size - in a sequence of sub- 
avalanches S — ^^ Sa- The number of sub- avalanches in 
the main avalanche is finite for a smooth landscape, and 
infinite for the continuous Brownian landscape. The to- 



tal size S 



- ui is however the same as for a = 0, due 



to the Middleton theorem. For instance in the landscape 
of figure [T]i, the main avalanche is divided into three 
large sub-avalanches, and for the continuous Brownian 
landscape the intermediate segments are also divided into 
smaller sub-avalanches, at infinitum. Figure [l] illustrates 
the correlation between the sub-avalanche structure (in 
u) and the realization of the random landscape, where 
larger hills favor the breakup into sub- avalanches. Note 
also that in intermediate regions where u is very small, 
rh starts decreasing again (it decreases whenever ii < h). 
The effective driving seen by the particle then becomes 
vn?w — aTh and increases. This mechanism triggers a new 
sub-avalanche, and so on. 

To obtain the dynamics one must solve the equations 



( 11 1, which we do below. For the standard ABBM model 
|46j . and in the mean-field theory of the elastic interface 
[T51 [5S] , it was seen that an avalanche terminates with 
probability 1, i.e. u{t) = for i > T. This allowed defin- 
ing and computing the distribution of avalanche dura- 
tions [IS SS] , and their average shape [551 ISSl US] • 

In presence of retardation, and for an exponential ker- 
nel, the avalanche duration defined in the same way be- 
comes infinite. Inside one avalanche, the velocity u{t) 
becomes zero infinitely often, but is then pushed for- 
ward again by the relaxation of the eddy-current pres- 
sure. Thus, an avalanche in the ABBM model with re- 
tardation splits into an infinite number of sub-avalanches, 
delimited by zeroes of u. Each sub-avalanche has a finite 
size Si and duration T^ with S — ^^Si (the same size as 
in the standard ABBM model), but ^ ■ Ti = oo. 

Below we study in detail two limits: 

In the slow-relaxation limit t ^ Tm the duration of 
the largest sub-avalanches remains of order r^, while the 
total duration is of order r. This leads to the estimate 
that the main avalanche breaks into ^ t/t„i significant 
(i.e. non-microscopic) sub- avalanches. 

For the fast-relaxation limit t <^ Tm {— ^ here) h fn ii. 
The correction to the domain-wall velocity u{u) is small 
in this limit, and vanishes as t/t^ — > (in contrast to the 
linfit r/Tm — > oo discussed above). In fact, the correction 



due to retardation amounts to a rescaling of the velocity 
as w — >■ (1 -|- aT)u. 

Of course, in presence of driving, the total duration is 
not strictly infinite since at some time-scale the driving 
will kick in again, and lead to another main avalanche, 
itself again divided in sub- avalanches and so on. We can 
call that scale again t„ but its precise value may differ 
from the estimate for the case a — Q. 

Thus one main property of the retarded ABBM model 
is that it leads to aftershocks, a feature not contained 
in the standard ABBM model. The main avalanche is 
divided into a series of aftershocks (the sub-avalanches) 
which can be unambiguously defined and attributed to 
a main avalanche (which basically contains all of them) 
in the limit of small driving. This sequence of sub- 
avalanches is also called an avalanche cluster. The after- 
shocks are triggered by the relaxation of the additional 
degree of freedom h. That in turn changes the force act- 
ing on the elastic system. Relaxation and aftershock clus- 
tering have been recognized as important ingredients of 
an effective description of earthquakes; the present model 
is a solvable case in this class. In some earthquake mod- 
els considered previously, relaxation was implemented in 
the disorder landscape itself [SUHSS] . Here the relaxation 
mechanism is simpler, which makes it amenable to an 
analytic treatment. Note of course that at this stage it 
is still rudimentary. First it is not clear how to identify 
the "main shock" among the sequence of sub-avalanches; 
while there is indeed some tendency, see e.g. Fig. [ijl, 
that the earliest sub-avalanche is the largest, this is not 
necessarily true. Second, to account for features such as 
the decay of activity in time as a power-law (Omori law 
[53]) one needs to go beyond the exponential kernel, to 
a power-law one. Finally, more ingredients are needed if 
one wants to account for other features of realistic earth- 
quakes, such as quasi-periodicity. 



B. Stationary motion 

In the case where the driving velocity is constant, 
w{t) = vt, the distributions of the domain-wall velocity 
ii and of the eddy current pressure h become stationary. 
The distribution of u for small ii has a power-law form 
with an exponent depending on v, 



Piu) 



-i+T 



(15) 



There is no contribution ^ 6{u). Vc is a critical driving 
velocity, which separates several different regimes: 

1. For V > Vc, the velocity ii never becomes zero. It is 
not possible to identify individual avalanches, one 
can say that there is a single infinite avalanche. 

2. For < V < Vc, the velocity ii vanishes infinitely 
often. The times {ti\u{ti) = 0} delimit individual 



(sub-) avalanches^ Their durations Ti :— ti+i — ti model, like for the quasi-static driving discussed above, 



and sizes Si = /^.'^^ dt'u{t') have distributions 
Pv{T) and Pv{S) d epend ing on the driving ve- 
locity V. In section [VIII we compute Pv{T) for 
the standard ABBM model and for a special case 
of the ABBM model with retardation. For sub- 
avalanches, starting at ui = 0, and a fixed value of 
the eddy-current pressure /ij, in the limit of small 
a and t — Tm, we show that 



PviT)^T 



-2+v+ahi 



for T ^ . 



In particular, the pure ABBM power-law exponent 
pS{T) - T-'^+-" is not modified for the first sub- 
avalanche, starting at hi = 0. Since the typical 
hi goes to zero as i; — >■ 0, we conjecture that the 
quasi-static exponents are still given by the mean- 
field values P{S) ~ S'^/"^, P{T) ~ T'^. 



In sections [V C| and |VI[ we compute Vc in several limiting 
regimes. For t ^ Tm, i-e. eddy-current relaxation slow 
with respect to the domain-wall motion, we obtain in 
section IV CI 



rilm^ + a) 



+ O [tJt) . 



This means that slow eddy-current relaxation decreases 
the critical velocity. The stronger the eddy-current pres- 
sure a, the smaller Wc becomes. On the other hand, for 
T <C T,„, i.e. fast eddy-current relaxation, we obtain in 
section I VI I 






Hence, fast eddy-current relaxation also decreases the 
critical velocity. However, the correction in this case is 
small and vanishes, as the time-scale separation between 
Tjn and T becomes stronger. 

The above regimes 1 and 2 do not change qualitatively 
compared to the standard ABBM model. This means 
that features like the power-law behavior of P{u) around 
M = are robust towards changes in the dynamics, as 
long as it remains monotonous. 



Non-stationary driving: Response to a finite 
l4icl4 



Instead of continuous driving, let us now perform a 
kick as defined in section II CI In the standard ABBM 



^ Note that there are no finite-time intervals where the velocity u is 
identically zero, since else the probability distribution | |15| l wrould 
have a S{u) part. Thus the times ti are single points, which may, 
however, be spaced arbitrarily close. Scaling arguments suggest 
that the set of points {U} has a fractal dimension of —. 



this leads to an avalanche on a time scale of order t„ 
which terminates with probability 1. At some time T, 
the domain-wall velocity ii becomes zero. The domain 
wall then stops completely, so that u{t) — for all t > T. 
This gives an unambiguous definition for the size and 
duration of the non-stationary avalanche following a kick 
|46| . Formally, this behavior is seen by computing the 
probability Pu(^t)=a- It turns out that Pu{t)=o > for any 
t > 0, and tends to 1 as t -^- oo. The distribution P{ut) 
(with Ut = u(t)) for t > has a continuous part and a 
(5-function part: P(ut) = Put=oS{ut) +'P{ut) [IS]- 

In the ABBM model with retardation, the situation 
is different. We show in section VD that p^, 



^Mt=0 



= 



following a kick, so that the dynamics never terminates 
completely. If one defines the avalanche duration T as 
T = ram{t\us — Ofors > t}, T is infinite. This is also 
seen from the example trajectories in figure [2b. How- 
ever, the velocity intermittently becomes zero an infinite 
number of times. Thus, the avalanche following a kick is 
split into an infinite number of sub-avalanches, just like 
a quasi-static avalanche discussed above. 

On the other hand, the sub-avalanches become smaller 
and smaller with time. In section VII A we show that 
the total avalanche size S := L dtiif following a kick 
of size Wo is finite and distributed according to the same 
law as in the standard ABBM model [45, 



P^oiS) 



Wo 



2\/TTaS2 



{wQ — m S) 
4ctS 



(16) 



This result holds independently of the memory kernel /. 
For infinitesimal kicks, wo — >■ 0, Pwa{S) becomes the dis- 
tribution of quasi-static avalanche sizes discussed above. 

The disorder- averaged velocity Ut following the kick de- 
cays smoothly. In the standard ABBM model, the decay 
is exponential |46j . With retardation, we show in section 
IV that the dependence of iit on t is directly related to 
the form of the memory kernel /. 

Another interesting observable is the mean avalanche 
shape. Conventionally, it is defined at stationary driving 
for a sub- avalanche: One takes two neighboring zeroes 
m(0) = and u{T) = which delimit a (sM6-)avalanche 
of duration T. The mean avalanche shape is then the 
average of the domain-wall velocity u{t) as a function 
of time, in the ensemble of all such (su6-) avalanches of 
duration T. It has been realized [22J that the skewness 
of this shape provides information on the relaxation of 
eddy currents. 

However, this definition is hard to treat analytically. 
Instead of considering the mean (sub-)avalanche shape 
at a constant duration, we discuss the mean shape of 
a complete avalanche (consisting of infinitely many sub- 
avalanches, with infinite total duration) of a fixed size S, 
triggered by a step in the force at i = 0. In section VII B 



we give an explicit expression for this shape at fixed size, 
for exponential eddy-current relaxation. We show how it 
reflects the time scale of eddy-current relaxation. 




2 4 6 

(bl) w{t) = 1.5, a = 0, T = 2, ri = 1, m = 1 



(b2) w{t) = 1.5, a = 1, T = 2, 77 = 1, m = 1. 




0.5 1.0 1.5 2.0 

(cl) w{t) = S{t), a = 0, r = 2, 77 = 1, m = 1. 



0.5 1.0 1.5 2.0 

(c2) w{t) = S{t), a = 1, r = 2, 77 = 1, m = 1. 



Figure 2. (Color online) Example trajectories for u{t) (thin, red) and h{t) (thick, yellow), for various parameter values. The 
left column (al),(bl),(cl) corresponds to the standard ABBM model {a — 0), the right column to the model with retardation 
(here a = 1). Figures (a), (b) correspond to stationary driving with a constant velocity, whereas the driving in (c) has a kick at 
t — 0. Observe that after a kick, ?i(i) in the standard ABBM model becomes zero permanently after a certain time, see figure 
(cl) , whereas in the ABBM model with retardation (c2) sub-avalanches restart infinitely often. 



The phenomenology discussed here is expected to be 
similar if instead of a kick at t = 0, one takes some arbi- 
trary driving wt for i < 0, which stops at i = so that 
Wt>o = 0. 

We see that the non-stationary relaxation properties 



of the retarded ABBM model differ qualitatively from 
those of the standard ABBM model. They provide a 
more sensitive way of distinguishing experimentally the 
effect of eddy currents than stationary observables at fi- 
nite velocity, and allow one to identify the form of the 



memory kernel /. In the following sections, we provide 
quantitative details underlying this picture. 



III. SOLUTION OF THE RETARDED ABBM 
MODEL 

In this section, we apply the methods developed in 
[T51 [5S1 Sni US] to obtain the following exact formula for 
the generating functional of domain-wall velocities. 



^/j At tit dt _ m /j WtUt dt 



(17) 



It is valid for an arbitrary monotonous driving wt > 0, 
where Uj is the solution of the following nonlocal instan- 
ton equation, 

/•OO 

r]dtu{t) - (to^ + a)u{t) + cru{tf ~a ds /'(s - t)u{s) 



^A(t), 



(18) 



with boundary condition u(oo) — 0. The important ob- 
servation that allows such an exact formula is that for 
monotonous driving, the motion in the ABBM model 
with retardation is still monotonous, as in the standard 
ABBM model (see appendix IA| as discussed above. 



To prove (17) we apply the same series of arguments 



as in the absence of retardation|THl [33 US]- Taking one 
derivative of Eq. ([5]) gives a closed equation of motion for 
u{t), instead oi u(t): 



ridtu{t) + au{t) + a / ds/'(t — s)u(s) 
= VW)^{t) + rr? \w{t) - iiify] . (19) 



^(t) is a Gaussian white noise, with ^{t)^{t') — 2a5{t — 
t'). The term ^u{t) comes from rewriting the position- 
dependent white noise in terms of a time- dependent white 
noise, 

2cr 



i{u{t))au(t')) = 2a5{u{t) - u{t')) = --5{t - t') 

u{t) 



au{t)) 



1 



m 



/u{t) 
dtF{u{t)) = um(u(t)) 



(20) 



This uses crucially the monotonicity of each trajectory. 

Using the Martin-Siggia-Rose method, we express the 
generating functional for solutions of ( 19 1 as a path in- 
tegral. 



o/t >^fiit — 



I?[u, u\e 



-S[u,u\+J^ Xtut 



S[u,u] 



Ut 



TjdtUt + aiit -\- a dsf'{t-s)u. 



■m^{wt - Ut) 



a I u^ut. (21) 



For compactness, we have noted time arguments via sub- 
scripts. We will use this notation from now on when 
convenient. 



As in the standard ABBM model, the action (21 ) is lin 



ear in ii. Thus, the path integral over ii can be evaluated 
exactly. It gives a (5-functional enforcing the instanton 
equation (18). The only term not involving ii in the ac- 
tion is rrrj utWt, which yields the result (17) for the 
generating functional. For more details, see section II in 
[Ml and sections II B-E in [18]. 



Similarly to the discussion in [T51I3B] the solution (17) 
generalizes to an elastic interface with d internal dimen- 
sions in a Brownian force landscape (i.e. elastically cou- 
pled ABBM models). There is indeed a simple way to 
introduce retardation in that model to satisfy the mono- 
tonicity property. We will not study this extension here. 

For the case of exponential relaxation, f{x) — e^^^'^, 
Eq. ( 19 ) can be simplified to a set of two local Langevin 



equations for the velocity u and the eddy-current pressure 
h: 



r\dtut = VUt6 + "1 {wt ~Ui\~a {ut - ht) (22) 

rdtht = ut- ht. (23) 

The action for this coupled system of equations is 



S[i 



■qdtUt + a(ut~ ht) + m"^ (ut - m) 



au?ut + ht (rdtht + ht - iit 



)} 



This action is linear in iit and ht- Thus, integrating over 
these fields gives 5-functionals enforcing a set of two local 
instanton equations for ut and ht, 

rjdtUt — (m^ -\- a) Ut + a-uj + ht ^ —^t, (24) 

rdtht ~ ht + aut = -A^t- (25) 

We then obtain the generating functional for the joint 
distribution of velocity ii and eddy-current pressure h, 



o/j(AtUt+Mt'»t)dt _ m^ f^wtutdt 



(26) 



in terms of the solution to these two instanton equations. 
It reduces to (17) for ^( = 0. 



Now, the remaining difficulty for arbitrary observables 
is to obtain sufficient information on the solutions of ( 24 ) 



and (25) with the corresponding source terms. We shall 



see that this is more difficult than in the standard ABBM 
model, but can be done for certain observables and cer- 
tain parameter values. 



A. Dimensions and scaling 

Before we proceed to compute observables, let us dis- 
cuss the scaling behaviour of our model, and determine 
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the number of free parameters. The mass m can be ehm- 
inated by dividing both sides of (191 by m^, 



■q 



dtu{t) 



Mt) 



dsdtfit - s)u{s) 



^W)(.it) + w{t) - Ht) 



(27) 



The time derivative -\dtu{t) shows that there is a nat- 
ural time scale r™ — r\lrr? so that t = t'T„i, where t' is 
dimensionless. The nonlinear term — 2\/u(t)£,{t) shows 
that there is a natural length scale Sm = -^, so that 
u = Smu' , where u' is dimensionless. We thus rescale 
velocities as 



u{t) 



ifm^ 



iu{t') , w{t) 



a 
rjm 



2^' it') 



(28) 



using the natural unit of velocity Vm — SraJTm- Multi- 
plying with vn?ri/a, we get the equation 

dru'{t') + ^ii'it') + ^ [ ds' dt'fit' - s')u'is') 



^^/i^e{t') + w\t')-u\t'), 



(29) 



where the noise is now (Citi)^' (h)) = 2S{ti — ^2)- Effec- 
tively, for the dynamics in terms of the primed variables 
we have m = a — i] — I (i.e. we have fixed the units of 
time and space so that Tm = Sm = !)■ 



For the standard ABBM model, a = 0, and Eq. ( 29 ) 



is a dimensionless equation without any free parame- 
ters. To describe a signal u{t) produced by the standard 
ABBM model, it thus suffices to fix the velocity (ampli- 
tude) scale Vm = ;j^: and the time scale t^ = :^- 

For the ABBM model with retardation, we have an 
additional time scale r, on which the memory kernel 
f{t — s) in (19) changes. The ratio of r to the time 
scale of domain-wall motion t„, = -^ is a dimensionless 



parameter r' := t JTm- Eq. (29 1 also contains a second di- 
mensionless parameter a! := -^j, which gives the strength 
of the eddy-current pressure, as compared to the driving 
w by the external magnetic field. We thus remain with 
two dimensionless parameters r' and a', which cannot be 
scaled away. 

From now on, we will use the rescaled (primed) vari- 
ables only. To simplify the notation, we drop all primes; 
we thus remain with the dimensionless equation of mo- 
tion 



dtiiit) + au{t) + a dsdtfit ~ s)u{s) (30) 

= VW)S.it) + w{t)-u{t). 

This amounts to setting m = a = ri = lm the original 
equation of motion, i.e. to working in the natural units 
for the ABBM model without retardation. 



IV. MEASURING THE MEMORY KERNEL / 

First, we discuss how the function / in equation ([s]) 
can be measured in an experiment or in a simulation. 
This allows verifying the validity of the exponential ap- 
proximation (l4|. We consider the mean velocity u{t) at 
t > following a kick by the driving field w{t) at i = 0, 
i.e. w{t) = wo6(t). Our claim is that its Fourier trans- 
form and the Fourier transform of the memory kernel / 
are related via 



die-'"* 



Wq 



Ht) = 2^. r" r, .y (31) 

TO^ + lui [i] + a J [lu)\ 



where /(w) := J^ dte-'"*/(i). 

To show this, we apply ( |17p to express the mean ve- 
locity at time to > a.s 



uito) ^ dx\^^^e^^(to) ^ dxl^^^C' 



J^dtuit)w{t) 



^d. 



A=0 



^wou{t=0;to) 



(32) 



The function u{t) is the solution of (18) with A(t) = 
XS{t — to). Since above we only need the term of order A, 
and u{t] to) is of order A itself, the nonlinear term in ( 18 ) 



can be neglected. In other words, the disorder does not 
influence the mean velocity u(io)i a-^d to obtain u(t; to)j 
it suffices to solve the linear equation 



r]dtu{t; to) - (m^ 
= -X6{t~to). 



a)u{t;to) — a ds f'{s ~ t)u(s;to) 

(33) 



Its solution is a function of the time difference t — to only, 
u{t; to) = u(t — to), which can be obtained by taking the 
Fourier transform u(uj) := J_ dr e~*"'^'u(r) as 



(iwri — TO — a)u{Lu) — au{uj) [~iujf{~uj) — 1] 



u{uj) 



A 



-iurj + TO^ — aiujf{—uj) 



-A 
(34) 



Here /(w) = L die *"*/(i) is the Fourier transform of 
the memory kernel. Inserting this relation into (32), the 



Fourier transform of the mean velocity after a kick is 



dio e-*"*''u(io) - Wo / dio e~^"*«u(0 - to) 



= Woui-Ld), 



(35) 



which then gives (31 1, as claimed. In fact, it is easy to 



see from ( 32 1 that a more general relation holds for a kick 

woiuj) 



of arbitrary shape. 



Ut 



m? -f iw [77 + aj{uj)] 



(36) 



where 100(1^) := /o°° die-'"*w(t). 

This relation allows one to obtain, at least in princi- 
ple, the memory kernel /(i) by measuring u[t) following 
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a kick. This permits to verify the validity of the expo- 
nential approximation ( 11 ) experimentally. It also allows 
to test the validity of the ABBM model. Indeed, while 



(36) at small wo(w) is simply a linear response, the fact 
that it holds for a kick of arbitrary amplitude is a very 
distinctive property of the ABBM model. Alternatively, 
it may allow to determine the frequency range in which 
the model provides a good description of the experiment. 



THE SLOW-RELAXATION LIMIT 



<r 



In order to go beyond the mean velocity and see the 
influence of disorder, one needs to solve the instanton 
equation (181 including the nonlinear term. Even in the 



special case of exponential relaxation, where ( 18 ) reduces 
to the local equations ( |24| ) and (25), their solution is 
complicated. However, we can analyze the latter in the 
slow-relaxation limit r^ — rj/ni? <^ t. In this limit, the 
relaxation of the domain wall to the next (zero force) 
metastable state, occurring on a time scale Tm, is much 
faster than the relaxation of eddy currents (occurring on 
a time scale r). Using the expressions for the relaxation 
times derived in (55], one sees that this is the case for 
very thick or very permeable samples^ To simplify the 
expressions, we rescale ii as discussed in section III A 
This amounts to setting m = a = rj — 1. Thus, the time 
scale of domain- wall motion becomes r ^ 1 . 

In the following sections, we will compute stationary 
distributions of the eddy-current pressure ht and domain- 
wall velocity Ut at constant driving Wt = vt, as well as 
their behaviour following a kick. A similar calculation 
for position differences at constant driving velocity is rel- 
egated to appendix IB] 



A. Stationary distribution of eddy-current pressure 



0), at constant driving 



Using ( 26 ) , the generating functional for the eddy- 
current pressure h = h{t 
Wt = vt is 



(37) 



u{t) is obtained from the instanton equations (24), (25) 
with the sources /i(i) — iJ,S{t), and A(i) — 0. From (25 1, 



one sees that h(t) evolves on a time scale s — t/r. On this 

scale, both h{t) and u{t) have a finite limit for -^ > 0. 

In this limit they are related via 

h{s) ^-u{sf + {l + a)u{s), 

u{s) ^Ua + 1- J{a + IY-Ah{s) ) . (39) 



(38) 



® In the pure ABBM model, the small-dissipation limit r; — > is 
equivalent — up to a choice of time scale - to the limit of quasi- 
static driving f — >■ 0"*". However, these two limits are different 
for the retarded ABBM model which we discuss here. 



The equation ( p5| for dsh{s) reads 

dsh(s) — h(s) — au{s) 



(40) 



Replacing h(s) on both sides of this equation using 
Eq. (38) yields a closed equation for m(s). 



[1 -I- a — 2u{s)]dsu{s) — u{s) — u{s)^ 



(41) 



The boundary condition at s = is fixed by the source, 
^(s) = /i(5(i) = 7(5(s) =: firS{s) (note u{s > 0) = h{s > 
0) = by causality): 

h{0) = Mr ^ "(0) = - (a + 1 - ^/Ja^TW^^^) ■ 

(42) 
Using Eq. (41 1, we can now compute the generating func- 

"(°) udu 



tional (37), 



u(t) dt 







u(s) ds 



dsu{u) 
"<°) {1 + a - 2u{s))du 
1-u 



= T [2w(0) + (1 - a) In (1 - u{0))] 
Inserting this result into Eq. ( [37| we get 



(43) 



(44) 



where mo(Mi) = ^(0) is given by ( |42[ ) and we have defined 
a rescaled velocity v^ '.— vt (i.e. the driving length during 
the relaxation time) 

The stationary distribution of h^ :— rh, obtained by 
inverting the Laplace transform, is 

f, [2vr~{?, + a)hrf\ 
X exp < hr - ^TT^ — — > X 



8/i 



X Ua- l)V2/lr£'(l-a)^,-l 

+2-0(1-0)1) 



(1 — a)hr + 2Vr 



(1 — a)/i,. + 2w,. 



(45) 



-D is a parabolic cylinder function [54 . 

For small /i^, the distribution (45) behaves as 



P{K) 



-=e «-r VrK - \ — 

v/TT V Vr 



{a-l)vr 



X (^l + 0{^r)) 



(46) 



Thus, there is a small-/i cutoff (due to the exponential 
term) and a power-law regime with a non-trivial expo- 
nent, hr ^ "^ ■ Note that the above results hold in 
the double limit w — >■ and r — >■ oo with v,, = vt fixed. 
Restoring units this is Vr = vt / Sm = t /t^ fixed, with 
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both T, Ty ^ Tm , hence v^ compares the two longest time 
scales, the driving time scale and the eddy-relaxation 
time scale. 

In the limit where the driving is slow compared to 
the eddy-relaxation time scale, Vr — > 0^, the station- 
ary distribution ( 45 ) takes the form of a limiting (un- 
normalized) density p{hr) proportional to 






1) < erf 



-(a- l)y/h 



1 



^1,3/2 



In this limit, the small-/ir behaviour is a pure power law, 
1 a- 1 



dvL^P{K) 



fnhr 



2hr 



0{hr 



(47) 



We note the resemblance of the tail of the distribution of 
h and the one of the size S in the usual ABBM model with 
the 3/2 exponent in both cases. If we assume that during 
avalanches (sub-avalanches) ii varies much faster than the 
relaxation time t (i.e. on scales r^ ^ r) we can rewrite 
'*(*) ~ Sat ^( S'Q.e~('~*°)/'^ where sub-avalanche a oc- 
curs at ta- Schematically h{t) integrates avalanche sizes 
occurring in a time window of order r, which could ac- 
count for the similarity. 



B. Eddy-current pressure following a kick 



Still in the limit -^ — > 0, let us now discuss a non- 
stationary situation: The dynamics following a kick of 
size Wo at i = to < 0, Wt — WQS{t — to). Using (26), the 
generating function for the eddy-current pressure at time 
is given by 



a = 


Us 


1 / 


= ^(l-\/l-4e''fio(l-ito)j 


a = 1 


Us 


= e^uo 


a = 3 


Us 




uie'l'^ - uoe"/ V^'oe"/' - 4{io + 4 


2(uo - 1) 



Table I. Some particular solutions of the implicit equation 
(|49[) describing the eddy-current pressure following a kick, in 



the limit ^- = 0. 

m T 



Taking the inverse Laplace transform, one obtains the 
distribution of eddy-current pressure h^ := t/i(0) after a 
kick of size wq at time to < 0, 



P{K) - 



Wo 



cto/(2-r) 



Wo 



2hr 



t^ 

_ (51) 

The average pressure hr = e*°/(^'^^wo/2 decays exponen- 
tially with time. Note that the limit to = 0~ leads to 
a non-trivial P(hr) which should hold within the entire 
matching region r^ ^ |to| ^ t. 

For a 7^ 1, we did not obtain an exact solution. How- 
ever, for any a, taking the limit /i — > — oo, or equiv- 
alently uo — > — oo, Eq. (491 shows that Ut — > — oo. 



This implies that the probability to find zero pressure, 
■ph=Q = limA^-oo e™""*!! — 0. Thus, after a kick at t = to, 
there is no time T > such that ht — for all t > T; the 
eddy current pressure never stops. A similar discussion 
for the domain- wall velocity follows in section [VD| 

Another interesting statement can be made regarding 
the time integral of the eddy-current pressure following 
a kick. From Eq. (fTTl) it must equal the total avalanche 
size (integrating this equation and using that for a kick 
h{0) = h{oo) = 0), i.e. J^ /i(t)dt = S. 



afJ,ho — g™OMto 



(48) 



Distribution of instantaneous velocities 



where Ut is the solution of (24 1, (25) with the sources 
A(t) = 0, fi{t) = fJ,S{t), as in the previous section. Now, 
we need its time-dependence and not just the total inte- 
gral. An implicit solution for u(t) at t < is obtained 
from Eq. (Ill]): 



du 



rdtu{u) 



■, , Ut , , , 1 — Ut 

= (1 + a) In — -h 1 - a) In -. 

uo 1 - Wo 



(49) 



Mo is fixed by (42). As in the previous section, we define 
a rescaled time s :— -. There is no expression in closed 
form for Ug for general a, but for specific values one ob- 
tains simple expressions (see table M. In the case a — 1, 
the solution is particularly simple. Equation (48) gives 



=M'i(0) 



oe'o/(2-)(i_yrr7i;:) 



(50) 



The distribution of instantaneous velocities P{u) at 
stationary driving is one of the simplest observables that 
can be determined from an experimental Barkhausen sig- 
nal. For the standard ABBM model, it has been obtained 
in |20) . For a d-dimensional elastic interface driven quasi- 
statically through short-range correlated disorder, this 
form is modified by universal corrections below the criti- 
cal dimension df. . These corrections have been computed 
using the functional renormalization group to one loop 
in e — dc — din [13 [5S]. Using (26), the generating func- 
tion of the instantaneous velocity, for constant driving 
Wt = vt, is 



(52) 
(|2l, 



Now Ut is the solution of the instanton equations 
( [25] ) with the sources A(t) = \S{t), ^i{t) = 0. 

To obtain the leading-order velocity distribution for 
T ^ 1, we need to solve (24) to order t~^. The solution 
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Ut, ht has two time scales, which are well-separated in 
the r — >■ cx) limit: i oc 1 and t (X t. We thus introduce 



s := t/T, 



and assume the scaling 



u{t)=:u^^^t) + T-^uf'\t)+OiT-^), \t\(xl (54) 
h{t)=:T-^U''^{t) + 0{T-^), |i|cxl (55) 



m 



:=■ -r-h~M) 



(,s) + 0(r-2), |i|(x 



^(t)=:T-i/i(")(s)+0(T-2), |i|oc 



(56) 

(57) 



Physically, the first regime |i| oc 1 is the regime where 
the eddy currents have not yet built up (h <C u) and 
are negligible. Hence the instanton is, up to a parameter 
change, identical to that of the standard ABBM model. 
The second regime is the regime of quasi-static relaxation 
of the eddy currents built up during the first stage. In 
that regime the instanton will be related to the instanton 
for the eddy-current relaxation discussed in the previous 
section. 

The source terms enforce the boundary conditions 
m('')(0) = -A, U'''>{0) = 0. We now construct u^'''''^ and 



1. Boundary layer: \t\ oc 1 

Let us first compute the leading term u^''\ which is of 
order 1. For -r < i < 0, inserting Eq. (l54| into Eq. ([25|, 



the term ht is subdominant compared to rdtht and Ut- 
We therefore obtain 



U''\t)^a / dt'u^''\t')+OiT-^) 



(58) 



Thus, the term h{t) in Eq. (24) is of order r ^, and neg- 



ligible in this regime. This is consistent with the inter- 
pretation of the boundary layer as the regime where the 
eddy currents have not yet built up. Eq. (24) reduces to 



dtu[''^-{l + a)u[''^+a{u[y = -Xt. 



(59) 



This is just the instanton equation (Eq. (13) of Ref. [5F) 
of the standard ABBM model, but with a modified mass, 
771^ = 1 ^ 1 + a. We obtain the known solution [551 SSI 



7iW(i) 



Consequently, for t 



(a+l)Ae("+i)* 
a + 1 + A (e(°+i)* - 1) ' 

-> — cx), ft,('')(i) is given by 



h^'''>{-oo) = a f dtu^''\t) = -a In (l —\ 



(60) 



(61) 



To compute the correction m]^ of order r ^, we need to 
expand (p4|) to the next order. We get the linear equation 



dtu'i\t) - (1 + a)uf\t) + 2u^^\t)uf\t) + /iW(f) = 0. 

(62) 



(53) n['>{t) 



Using the expressions (58), (60) for /i'^^'(t), its solution is 
given by 



(l + a)[l + a + (l-e-(i+'^)*)A]' 



- Ae^'^+i)* 



X Li 



2(l + a- A)x 



Lij 



1+a-A/ V 1+a-A 

- (1 + a)t{l + a - A) + A (e^'^+i'* - l) 
(l + a-A)2-A2e2(-+i)*" x 

xln 1 + ^ '- 

+2{a + l)Ate('^+i)*(l + a - A) In (l —\ \ 



For t — > — oo, this tends to a constant, 

^-oo 1 ^ ^ l + a \ l + aj 






(63) 



(64) 



Since u^ '{—oo) = 0, see Eq. (60), this is the dominant 



contribution of the boundary-layer solution for u in the 
limit t -^ — cx). 



2. Long-time regime: \t\ oc r 

Now, let us consider the regime t < — r. Inserting the 
rescaled time s := t/r into Eq. (24), we see that the 
term dtu{t) = T^^dgU^s) is subdominant in t~^. The 
instanton in this regime is thus a special case of the in- 
stanton discussed in section V A Applying Eq. ( 56 ) we 



see that u{t) oa t~^ is small. Thus, we can also neglect 
the non-linear term u{t)'^ in (24). This gives the simple 
relation 



l + a 



{(W^=^- 



Consequently Eq. (25) reduces to 



l + a 

The boundary condition at s = is now non-trivial, and 
given not by the sources, but by the asymptotics of the 
boundary layer as t — > — cx). 



h^^^O) = /i(^)(-oo) = -a In (l —\ 



The resulting solution of Eq. ( 24 ) is 



/i(°)(s) = /j(°)(0)eTTs = -a In (l —\ 



gl + a 



«(")(«)=- 



1 -l-a 



;T+^ In 1 - 



A 



a + 1 



(65) 

(66) 
(67) 
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Figure 3. (Color online). Instanton solution ut, ht of (241, 
( 25 I with sources \{t) — —S{t), /i(f) = in the slow-relaxation 
limit r 2> 1. Parameters are a — 1.3, r = 2. Yellow (thick) 
curve: ut, red (thin) curve: jr^ht, black dotted curve: long- 
time asymptotics (661, grey dashed curve: short-time asymp- 



totics ( |60[ ). The inset shows details of the boundary layer 
Itl oc 1. 



A non-trivial consistency check is that this expression 
matches the 0{t~-^) term of the i — )■ — oo asymptotics of 
the boundary layer given in Eq. (64), 



uW(0)=uf)(~oo). 



(68) 



The boundary-layer solution (roO| for t ex 1, and the long- 
time asymptotics ([66]) for i ex r compare well to a di- 
rect numerical solution of ( 24 ) , ( 25 ) in the corresponding 
regimes, see figure [3j 



3. The velocity distribution 

From the combined knowledge of the previous sections, 
we can extract the generating function for the velocity 
distribution ( 52 ) . We have 



Ut 



In 1 



X 



a+l 



aln 1 



a + l 



(1 + a) In 1 



A 



a + l 



(69) 



Thus, the generating function for the velocity distribu- 
tion (52) is 



sAtt — 



1- 



A 



a+l 



-{l+a)v 



i + o{Ty 



(70) 



The Laplace inversion is easy to do, giving the distribu- 
tion of instantaneous velocities to leading order in r^^ 
(but without any approximation in v). Restoring units. 



this is 



Piii) 



u— [m + a ) 



r[2(m2 + a)w] ii 



(71) 

We can compare this to the distribution of the standard 
ABBM model, 



P(u) = 



V a 



FT" i"~ 



m 



(72) 



We see that the effect of eddy currents on the instanta- 
neous velocity distribution, in the limit of r -> oo, is the 
same as if the mass in the standard ABBM model were 
increased from m^ to m'^ + a. In particular, this means 
that the transition between intermittent avalanches and 
continuous motion happens for driving velocities v re- 
duced by a factor of (1 + a/ni^). In dimensionful units 



(73) 



ri(m^ + a) 



D. Velocity following a kick 

Let us now assume that the driving velocity under- 
goes a kick at time i < 0, i.e. w{t') = WQ6{t' — t). As 
discussed in section I C we consider an initial condition 
at i = prepared in the "Middleton attractor" u{w) 
with ■a(O) — h{0) — 0. The kick gives deterministically 
u{t^) = rri^WQ =: Mj, so that the distribution of velocities 
iif := iiQ is the propagator V{u{,0\ui,t) at zero driving 
velocitjF] Applying (17), the generating function for the 
velocity uq at time t = is given by 



(74) 



One can obtain the probability distribution of uq in the 
slow-relaxation limit r 3> 1 by inverse Laplace transfor- 
mation as follows. There are two time regimes: 



(i) for \t\ ^ 1 we need to insert Ut 



nW 



in (74) as given by (60) and (63). For small t we need 



w- 



-i,~,W 



(t) 



to use only ( 60 ) and we recover the velocity distribu- 



tion and propagator of the standard ABBM model at 
zero driving velocity, but with modified parameters. The 
propagator of the standard ABBM model has been dis- 
cussed, among others, in [24l |25l |46]. Here we use the 
result from Eq. (24) in Ref. [25], or equivalently (19) in 



"^ In this section, we still keep r) 
some places for clarity. 



1 but restore units of m in 
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Ref. [IB], with V — and rn^ 
P(u) =P(u;0|m^wo;t) = 



'2sinh[it(a + m2)] 

\ g(a+m2)i 



X exp ■ 



X/l 



sinh (i (to^ + a) f) 



(75) 



This distribution is not normaHzed; formally, there is an 
additional (5- function term as in Eq. (24) in Ref. [25] . 
This indicates that this formula is valid for w ^ 1 only; 
there is another regime ii ^ 1/t <C 1 which requires a 
more careful treatment but will not be considered here. 

At later times there is a complicated crossover to 
regime (ii), which requires to keep the 1/r correction, 
which we will not detail here. 

(ii) for long times |i| ~ t, ut is given by (56) and (67 1. 
This gives 



aAllo 



1 



A 



1 



'""0 exnl" ^ "1 



(76) 



Inverting the Laplace transform, one obtains the velocity 
distribution and propagator 



P{u) 



1) 



ctp-it(l+a) 



T{et)u'^ 



£t 



awo 



l)r 



e-("+i) 



(77) 

It approaches rather quickly a (normalized and regular- 
ized) power-law distribution proportional to l/ii. Note 
that for an infinitesimal kick we recover the limiting 
density obtained from the stationary motion p{u) ~ 



where P^*^^ (ii) was obtained in ( 71 1 and rtyp is a typ- 
ical time scale. 

In general, the velocity distribution P{u) following a 
kick can be decomposed as 



Piu) = Pu=(i5{u) + Vrcgiu) , 



(78) 



where Pu^o is the probability that the domain-wall has 
come to a complete halt. From (77) one sees that in 
the ABBM model with retardation, Pu=o = and the 
domain-wall motion following a kick never stops com- 
pletely. This is in contrast to the standard ABBM model, 
where one has (cf. gB], Eq. (28)) 



Pu=o = exp 



Wo 



-1 



This can also be seen directly from the instanton solution 
u: The decomposition ( 78 ) implies 



oXii — 



e^"P(u) 



A— > — c 



-> Pu=0 



(79) 



Since e-*'" = Q^out^ ^^ ^g^j-^ conclude that Pu=o is zero, 
if and only if \un.\^^c>oUt = — oo. This is the case 



in the retarded ABBM model, where (67 1 shows that 

A— f — oo 



Ut oc 



-In 



('-i^) 



-> — oo. However, it is not the 

oUt ^ 



case in the standard ABBM model, where limA-j 
j^ is finite (cf. [46j, Eq. (14)). 

We conclude that in the ABBM model with retarda- 
tion, the velocity following a kick never becomes zero 
permanently, even though its mean decays exponentially 
in time over time scales of order r, with for t < — r. 



{uo) = 



awo 

(a + l)2i 



(80) 



Although the calculation above was done at leading or- 
der in T~^, we expect the phenomenology to be similar 
for arbitrary r. To make this explicit, we now consider 



the opposite limit of fast relaxation, 
analytical progress is also possible. 



T < 



in which 



VI. THE FAST-RELAXATION LIMIT r < ^ 

We can also consider the limit r ^ t„i = i]/m^, where 
the eddy currents relax much faster than the domain- wall 
motion. Experimentally, this limit is even more relevant 
than the slow-relaxation limit discussed in section |Vl As 
a function of sample thickness 6, the eddy-current relax- 
ation time T (X b^, whereas the domain- wall motion oc- 
curs on a time scale -^ oc ab. For typical experimental 
setups [22l |55] a:^ b and hence r <C ry/m^ . 

We now discuss the stationary velocity distribution, 
and the velocity following a kick in the driving velocity, 
in the fast-relaxation limit. As in section [V C[ we need 
to construct the instanton u, h solving Eqs. (24 1, (25) 
with sources A(t) = XS{t), fi{t) = 0. Now, however, r 
and not t~^ is a small parameter. We expect a two-scale 
solution: A boundary layer for \t\ ex t around t = 0, and 
an asymptotic regime for |t| ex 1. We thus introduce the 
rescaled time s := t/r and make the ansatz 



u(i) 
hit) 
u{t) 

m 



ui'^ it) + Tui'\t) + 0{t)\ |i|cxl 
h^i\t)+0{r), |i|cxl 

u\^\s) + TU^^\s) + 0{Tf , |i|(XT 

U-^\s) + 0(t), |i|(xr. 
A. Leading order 



(81) 
(82) 
(83) 
(84) 



At order r°, the instanton equations (25) and (24) re- 
duce in the asymptotic regime to 



dtii'i'^ (t) - (1 + a)u^^^ {t) -f \u^^^ (i)l + h^^^ (t) = 



-Ab) 



hl,"'{t) + au'o>{t)=0 



r.ib). 
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with boundary condition uo(0) = A. This gives the 
leading-order solution 






(85) 



In the boundary layer, the corresponding solution is 
u^o\s) = u^o\0) = A, and Eq. pll) gives 



;(a) 



(a). 



d,hi;'>{s)-h}j''>is) + aX^O 



B. Next-to-leading order 



We obtained in Eq. ([85| the leading-order solution 
Uo(t), valid in both regimes. Expanding around it, set- 
ting u{t) = uo{t) + TUi{t) + 0(r)^, we get an equation 
for ui(t) 

dtUi{t)~{l+a)ui{t)+2iloit)ui{t) + - \h{t) - auo{t)\ =0. 

T L 

In the boundary layer, dtUi{t) = -dgUi (s), and 

/io(s) — auo(s) = aA (1 — e*) — aX — —aXe". 

Hence, the next-to-leading-order contribution ui' (s) in 
the boundary layer satisfies 



r,W 



dsu\''>{s)=aXe' 



„--,W 



ui"^(s) = -aA(l-e") 



(86) 



On the other hand, Eq. (251 gives in the asymptotic 
regime 



Inserting this relation into Eq. p4| ) gives 

dtuf\t) - nf\t) + 2u'^^\s)uf\t) + adtu^*'\t) = 

-(b)., aXe\Xt-t-l) 
^ "i ^^' ~ (Ae* - A + 1)2 ■ 

Here we used the matching condition u\ ' (0) 
ul" (—oo) = —aX, as given by Eq. (86). 



(87) 



The next-to- leading order corrections ( 86 1 and ( 88 ) 



compare well to a direct numerical solution of Eqs. (|24|) 
and (25 1, see figure |4J 



C. Stationary velocity distribution 

With the above analysis, we can obtain some results on 
the velocity distribution. The integral over the instanton 
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Figure 4. (Color online) Instanton solution u, h of Eqs. ( |24[ |, 
(251 with sources X{t) — —5{t), jj,{t) = for fast eddy-current 
relaxation, r <C rj/rn^. Parameters are o = 5, r = 0.1. Yel- 
low (thick) curve: u{t), red (thin) curve: -h{t). Black dotted 
curve: Leading-order result (85 I, corresponding to the stan- 
dard ABBM model. Grey dashed curve: Solution (851 plus 



next-to-leading order correction ( 88 1 in the long-time regime 



|f I oc 1. The inset shows details of the boundary layer \t\ oc r. 
Blue dot-dashed curve (inset only): Solution ( |85[ ) plus next- 
to-leading order correction ( 86 1 in the boundary layer. 



solution gives 



u{t)dt = / dt u'-o^ {t)+T dt uf^ (t) + ©(r) 



ln(l - A) + ar 



X-1 



ln(l - A) 



(89) 



Inserting this result into Eq. ( 52 1 for the generating func- 



tional of instantaneous velocities gives 

^ - (1 - A)""('+°"^ ga.r^+0(r)^_ (gg) 

To the same order in r, this can also be rewritten as 

A 



u{t)dt = - (l + aT)ln 



1- 



^Xu — 



1 - 



A 



(1 + ar) 



l + ar 

-v{l+aT) 



+ 0(r)^ (91) 



0{t)' 



This makes it clearer that, to leading order, the form of 
the velocity distribution is not modified, and only the 
parameters are rescaled. 

For small u, it indicates that, in the fast-relaxation 
limit, the instantaneous velocity distribution P{u) has a 
power-law behaviour 

P(u)^u^i+"(i+-+'^M'). (92) 

Putting back the units, the power-law exponent becomes 
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-1 + — , where 



rjm-^ 



1-0-+© —-^ 



O 



'm? [rj -\- ar) \ri/m 



(93) 



We see that fast eddy-current relaxation decreases v^, 
just as in Eq. (71 1 for slow eddy-current relaxation. Both 
formulas have the form Vc = ^21 "^ar *) ' '^here Tf is the 
fastest time scale in the problem, Tf = rj/m? for the slow- 
relaxation limit, and Tf = r for the fast-relaxation limit. 



In contrast to Eq. (71 1 however, the correction we obtain 
here is perturbative: It vanishes as r — )■ 0. In the limit 
T -> we recover the standard ABBM model. 



D. Velocity following a kick 




Figure 5. (Color online) Instanton u{t) at a fixed time t = —1, 
as a function of A, for r = 0.1 and a = 1. Thick yellow line: 
Numerical solution of ( |24[) , (25 1. Thin red line: Leading- 
order, ABBM solution ( |85[ l, corresponding to a = 0. Dot- 
dashed black line: Next-to-leading order correction (88 1. 



The generating functional for the distribution of veloc- 
ities m(0) following a kick of size w dX t < Q can also be 
expressed in terms of the instanton solution (81 1, 



^XibQ 



e»"(t) ^ exp 



~(b)/ 
WUq ' 



(t) + wTuf\t) + Oir)' 



.~.ib) 



(94) 
and u^' ^^^6 given by Eqs. (85 1 and (88) above. They 



^-.W 



have a finite limit as A 



Hm u^'') = J 

A-y-00 1 — e ' 



A^-c 



uf^ 



ate 



(l-e-*)2 ■ 



(95) 
(96) 



This would suggest, that the velocity distribution con- 
tains a term ~ 5{u). However, for large negative A, the 
expansion above breaks down, and higher orders in r be- 
come non- negligible. By solving the complete instanton 
equations numerically one obtains figure [5J One observes 
that the leading order (standard ABBM) instanton ( 85 ) 
goes to a fixed value for A — >■ —00. The next-to-leading 
order correction ( 88 ) coincides better with the numerical 
solution, but still breaks down around A « —10, and goes 
to a fixed value, too. However, the true (numerically ob- 
tained) solution of the instanton equations goes to —00 as 
A -^ -00. Hence, limA->-oo e^"o = limA^-00 e"'"^*) = 0, 
and the distribution P{ua) does not have a S{ua) piece, 
consistent with the results obtained above in section FVDI 
in the r — )■ 00 limit. 

From the instanton expansion ( 85 ) , ( 88 1 , valid for A oc 
1, we can obtain the velocity distribution P{uo) following 



a kick at i < 0, in the regime up ct 1- Using Eq. (94), we 
write its generating function to order t as 



oXuq 



= cxp wA — 



•iC 



wB 



\-q {\-qy 



(97) 



with 



1 

1-e* 



A=- 



ae t 



,t 



1 



(e* - 1)2 



^=(^^31)^+^- 



[e\t + l)+t-l] 
(e* - 1)3 



C^T 



ae*[e\t + l) - 1] 
(e* - 1)4 ■ 



The inverse Laplace transform of e^"" can be written as 



Piu) 



wA—qii 



(98) 



27r 



-^ re**^ exp (w 

ZTT \ 



rCe-2# ^^-i4, 



-ure 



We have set X = q + re*"^ to arrive at the above formula. 
We numerically checked that the integral is independent 
of r. One can evaluate it analytically, if either B = 0, or 
C = 0, by expanding in powers of r, and retaining only 
terms which scale as r° (note r" ~ e'*^"). The final result 
can be written as a convolution of the two: 



P{u) =e 



wA — qu 



duiPi{ui)P2{u~ fii). 







Piiii) =d{u) + \l -T^/i {2^/Bv 



Bw 



P2iu)^Siu) + CuwoF2[ ^,2 



-Cu^w 



(99) 

(100) 
(101) 



Note again that formally the (5-function parts arc an ar- 
tifact of our expansion, which is not valid as A — )• —00 
or -u — >■ 0. We expect them to be smeared out on a scale 
e*/'^, which goes to as i — >■ —00. However, physically. 
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Figure 6. (Color online) Velocity distribution P{uo) given 
by (991 (without the 5-part), after a kick at time t < 0. 
The parameters were chosen tobew = l,a = l,r = 
0.2. The time of the kick varies from right to left as i = 
-0.3,-0.4,-0.7,-1.5,-2.5. 



these velocities are extremely small and unlikely to be 
observable. Thus the ^-function term is physically sen- 
sible, and can be interpreted as the probability that all 
significant avalanche activity has stopped. 

Numerically, the convolution can easily be computed. 
An example of the distributions for various times is shown 
in figure [6) We see that for small times the distribu- 
tion is peaked around the value w imposed by the step 
in the force. Later on the typical value of the velocity 
approaches 0, and the distribution becomes monotonous. 
Its area decreases since part of the probability is absorbed 
by the (smoothened) ^-function near m = (A = — oo), 
which we are unable to analyze here in more detail. 



VII. AVALANCHE STATISTICS AT FIXED SIZE 

In the previous section we saw, at least in the two lim- 
its 77/m^ <^ T and rj/m^ ^ r, that avalanches following 
even an infinitesimal kick never completely stop. Com- 
puting observables conditioned to their duration of first 
return to u = 0, i.e. the sub-avalanche duration, requires 
introducing an artificial "absorbing boundary" at w = 
which will terminate the avalanche once ii becomes zerco 
This task is deferred to section [VIII[ However, the mean 
velocity following a (finite or infinitesimal) kick still de- 
creases, and the total avalanche size remains finite. We 
will now compute its distribution, and other observables 
conditioned on the total avalanche size. 



Avalanche sizes 



We define the size S* of a non-stationary avalanche fol- 
lowing a kick of size wq at t = as S = J_ u{t)dt. The 
Laplace transform of the probability distribution of S is 



given by Eq. (17), 



aXS - p^/~oo'i(*)dt 



(102) 



Here Ut is the solution of ( |18[ ) with a time-independent 
source X(t) = A. This means that ut — u is also time- 
independent. Then, using /(O) = 1 and /(oo) ~ 0, the 
terms proportional to a drop out from Eq. ( |18[ ) and we 
get 



m u + au 



-A. 



Choosing the solution which tends to as A 

™2 



^r, 



4Acr 



and 



2ct 



2(T 



Inverting the Laplace transform gives 



P{S) = 



Wq 



iaS 



(103) 
0, we get 

(104) 
(105) 
(106) 



with S = Wq. Note that this extends to any finite kick 
of arbitrary shape replacing wq — J^ dtw{t) [TH]. This 
is precisely the distribution obtained for the standard 
ABBM model and the mean-field theory of interfaces in 
]IE[ 1351 H^. Of course, this can already be seen from 
the fact that the terms proportional to a drop out from 
(fT8| when u is time-independent. Note that this result 
is independent of the shape of the memory kernel / in 
([5]). This is a consequence of the monotonicity of the 
model, as discussed in the Introduction. In the limit 
of an infinitesimal kick, i.e. small wq, one recovers the 
stationary avalanche-size density. 

Universal corrections to the distribution ( |106[ ) are ex- 
pected when one goes beyond the mean-field limit and 
considers d-dimensional elastic interfaces. Without re- 
tardation effects, the universal corrections at slow driv- 
ing were obtained to one loop in an expansion around 
the critical dimension [TH]. We expect them to remain 
unchanged by retardation effects, as seen in this section 
for the mean-field case. 



B. Avalanche shape at fixed size 



The natural boundary at -ii = would be reflecting, since if ut 
becomes zero at some instant of time, it immediately restarts 
to positive velocities due to the decrease of the eddy current 
pressure in the next time step. 



The avalanche shape is usually obtained by computing 
the mean velocity as a function of time, in the ensemble 
of all avalanches of a fixed duration [1^1 HH [Ml I5S] ■ Here 
we shall instead consider the ensemble of all avalanches of 
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a fixed size S. In a numerical simulation or in an experi- 
ment, the shape at fixed size is just as easily measurable 
as the shape at fixed duration. However, it is easier to 
obtain theoretically with our methods, and it can be de- 
fined without a microscopic cutoff. We will thus compute 
the shape function defined via 



sit,S):-- 



diit iit Piut\S) 



J^ dutUtP{ut,S) 
PiS) 



(107) 
P{S) is the avalanche size distribution ( |106[ ). 

We follow the approach used in [T5J 1451 146] to obtain 
the avalanche shape in the standard ABBM model from 
the Martin-Siggia-Rose field theory. The driving wt per- 
forms a kick at i = 0, i.e. we set wt — woS{t). Wc then 
consider the observable 



s(to, A) := ut„ e^s = a^|^^^e^^+^'"o 



(108) 



s is related to the shape function s via a Laplace trans- 
form, 



s(io,A) 



dSs{t,S)P{S)e 



\s 



5 as defined in (108) can be evaluated using (17) 



s(to,A) 



^^1^=0^" 



(109) 



(110) 



where uo{fi) is the solution of (181 with the source At = 
X + fiS{t — to). To compute (110 1, we need to solve the 



instanton equation ( 18 1 to first order in ii. The solution 



for /i = is the constant u{X), obtained previously in Eq. 
( 104 1 for the size distribution. The correction of order /i, 
u'^^' , has to satisfy the linear (but still non-local) equation 



a ~(i) 



dsf{s-t)ui'^^ 

= -fiS{t-to). (Ill) 

We now restrict ourselves to the case of an exponentially 
decaying memory term, f{t) = e"*'"^. Through the sub- 



„"-,(!) 



stitution ul = ji{l — Tdt)gt, equation (111) is trans- 
formed into a linear second-order ODE, 

{dt-l-a + 2u){l - Tdt)gt + agt = -S{t - to). (112) 

The right-hand side yields the boundary conditions 
g{to) = 0, g'{to) = — 1/t. The resulting solution for 
gt is 



2re^t-to)/r 



9t 



2a(t-to) 



r{t-to) 



Aar ' 



The shape function (110) is then 

s(io,A) = e'^'"'"(^=")wou^^^(0) 



= e'^fi-^i^^^wo [.9(0) - Tg'(O)] 



^^(l_yT34A)^^ 



4ae" 



Aar + r^ 



i!^z!!^e-*«/^ 



(114) 



The shape at fixed size S is finally obtained by invert- 
ing the Laplace transform. This is best done using the 
coordinate r introduced in Eq. (113): 



dutUtP{ut,S) 



tq H-ioo 



dr 

27ri 



it,X{r))c 



'A(r)S 



dA 
dr 



D-t-ioo j^ (r- 2aT)(r + 2) r 



Aare r -f re 2t 



27ri ^r'^r^ 

S'(r + 2)2(r-2aT)2 _S _t 

(2 + r)(r-2ar) 



X cxp 
X Wo exp 



Wo 



1 



2rr 



(115) 



To fixes the location of the integration contour; it can be 
chosen arbitrarily, as long as r-g 7^ 1^ 

Our final result for the avalanche shape at fixed size 
(following a kick of arbitrary strength wo) is thus ob- 
tained by inserting this result into 

s(t,S) = 2y/^S^''^e 4S — / dutUtP{ut,S), 

Wo Jo 



where we have used (106 1. On this expression the limit 



of Wo —^ is easy to take and provides the result for the 
stationary avalanches. 

One non-trivial check of this formula is that the result- 
ing shape is properly normalized, 



dtsit,S) 
Jo 



' "* P(S) 



dr (r^ -I- 4ar) 

S'(r + 2)2(r-2aT)2 
IGr^r^ 



X exp 



3 (wQ-sy 
X S" 2 e 4S exp 



Wo {2 + r){r-2aT) . 

2 ^ 2rT ' 






Ai^/TT 



where in the second line we used the substitution p 

{r+2){r-2aT) 



where we defined r via 

(2 + r)(r 



2ar) 



2rr 



V1-4A. 



(113) 



^ This is in order to avoid the singularity of the integrand at r = 0. 
We checked that the result of the (numerical) integration of | |115[ l 
is independent of the value and the sign of rg . 
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<ii(t)\S> 




Figure 7. (Color online). Avalanche shape at fixed size 
S — 0.8, T = 0.5 and wq = 0. Curves are, from top (black 
solid line) to bottom (blue dot-dashed line), a = (standard 
ABBM model, as given by Eq. (|116[)), a = 1, a = 2, a = 5. 



The integral (1151 could be calculated in closed form 
in the case a = of the standard ABBM model. There 
one finds 



dutUtPiut,S) = 



(2t + wo)woe' 



S^^2S^0 + {2t + ^0)^ 



27^53/2 



and the result for the shape, 

s(i,5) = (2i + wo)e"*'*+'"°^/^. 



(116) 



In the limit of an infinitesimal kick uiq — > we thus obtain 
the shape at fixed size for the standard ABBM model 
(a ~ 0) for stationary avalanches as 



s(i,5) = 2te-*'/^. 



(117) 



For a > we could not find a closed expression, however 
the integral (1151 is easily evaluated numerically. Some 
example curves are shown in figure [7] in the limit of small 
wq. Observe that especially for large values of a, the ad- 
ditional time scale introduced by the eddy-current relax- 
ation is clearly visible. Overall the shape stretches longer 
in time, and becomes non-monotonous, as a is increased. 



1. Tail of the shape function 



The behaviour of the avalanche shape for long times 



can also be understood analytically from equation ( 115 1. 
For simplicity, we consider the case ujq = in the follow- 
ing. For large t, the integral is dominated by its saddle- 
point. Since we have e^^ 3> e^^r for all times, the 
dominant exponential factor is 



e«('-) := exp 



5'(r-f2)2(r-2aT)2 S t 



2a A 



Its maximum for large t is obtained by solving H'{r) = 0: 



2atT^ 



S 



+ -{aT-l) + 0{t) 



-1/3 



Determining the location of the saddle-point to higher 
order in t is more complicated. The terms of ©(i)"^/^ 
depend on whether the sub-exponential terms in ( |115 ) 
are included in the maximization procedure or not. How- 
ever, to the order given here, r,„ is independent of such 
choices. 



Since the integral (115) does not depend on rg as dis- 
cussed above, we can choose vq = rm- Setting r = u + iv, 
we have d^H{r) = —d^H{r) due to the Cauchy-Riemann 



equations. Thus, we can approximate the integral (115) 
for large t and fixed S by 



s{t,S) 




27r (r„ - 2ar)(r™ + 2)a e^(''-) 
Anrrl P{S) 

t 3 fatV'^ fS^ ^'^ 



T 2 V T 



at\ 1 — ar f S\ 



T I T \2j 

S (a^T^ - 5aT + 1) 



6r2 
5^/3a2/3 

V3(2i)l/3r2/3 



2 

2/3 



0{t) 



-1/3 



o{t)-^^^ 



(118) 



Again, note that to this order both the exponent and the 
prefactor are independent of whether the sub-exponential 
terms are included in the maximization. In particular, 
the term of order 0(i)° in the exponent is independent 
of the term of order 0{t)~^^^ in rm- We thus see that 
the Gaussian tail of the shape in the standard ABBM 
model is replaced by an exponential tail, decaying on a 
time scale r. This is confirmed by numerical Laplace in- 
version of Eq. (115), see figure [8J We also observe good 
agreement between the asymptotic expansion and the nu- 
merical result. 

Using a similar method one can try to determine the 
tail of s{t, S) for fixed t, at large S. In this limit, the 
maximum obtained by solving H' (r) = is 



2aT + 0{Sy 



One finds 



H"{r„ 



_S 

(1 + ar 



o{sy 

2 



-s + o{sf. 



Noting that P{S) ^ e~^'^, this means that the exponen- 
tial factor in the saddle-point contribution to 5{t,S) ^ 
e^'-'^'^') /P[S) vanishes to leading order. This indicates 
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equations (|T7|) and (26) for the generating functional of 



0.01 ■ 



0.001 ■ 



Figure 8. (Color online) Tail of the avalanche shape at fixed 
size S = 0.8. Blue dotted fine: Gaussian ABBM tail ( [Tl6 l. 
Red (thin) and yellow (thick) lines: numerical Laplace inver- 
sion of (1151, for T = 0.7 and r — 0.4, respectively. Black 
dot-dashed and dashed lines: Asymptotics (1181 for the cor- 
responding values, a = 1 in all cases. 



that s{t, S) will scale as a power-law for fixed t at large 
S. However, since the pre-exponential factors in (115) 
also vanish at r™ = 2ar, obtaining a quantitative result 
requires a more controlled approximation. 



VIII. SUB-AVALANCHE STATISTICS 



As we saw in section [VD[ an avalanche in the ABBM 
model with exponential retardation never strictly termi- 
nates, even after the driving has stopped. It is thus in- 
teresting to explore the "sub"-avalanches, or aftershocks, 
inside an avalanche, and their durations Ti, z = 1, .... Ti 
is defined as the time it takes u to start from u = at 
time ti^i = '^j=iTj, go to positive values u > and 

back to M = at time ti — X^i^i-T;'' without touching 
M = in the interval ]Oti_i,ti[. In other words, it is the 
separation in time between successive passages of ii at 
zero. The same question can be asked for avalanches at 
non-zero driving velocity ?; > in the standard ABBM 
model, which can also be seen as sub-avalanches of an 
infinite avalanche (there too u never vanishes on a finite 
time interval) p°] We will obtain detailed results in that 
case. 

A convenient setting to study this problem is the 
Fokker-Planck approach, introducing an artificial ab- 
sorbing boundary at u = 0. It can be implemented in 
the case of the simple exponential relaxation (|8| which re- 
duces to two coupled Langevin equations for ii and h ( 11 ). 



We will discover that, with such an absorbing boundary, 



^^ In the context of the standard ABBM models, 
avalanches are also called pulses |57| . 



these sub- 



domain-wall velocities, as well as some details of the in- 
stanton method, need to be modified. 



A. Sub-avalanches in the standard ABBM model 
at finite driving velocity 

In order to present our approach on a simple exam- 
ple, let us consider first the standard ABBM model with 
monotonous, but otherwise arbitrary driving w(t) > 0. 
The equation of motion ( 19 ) with a = is, due to its 



Markovian nature, conrpletely characterized by the prop- 
agator 



Villi; tf\ui;ti) := E[S{utj - ut)\ut, = Uj] 



(119) 



Here E means "expectation", i.e. the average over the 
disorder. As a function of the final velocity iii, Pt{u) — 
V{u[ — u;t[ = t\ui;ti) satisfies the forward Fokker- 
Planck equation 

dtPt{u) = dluPtiii) + du{u - Wt)Pt{u) = duJt{u), 

(120) 



where 



Jt{u) := duuPt{u) + (it - Wt)Pt{u) 



(121) 



is (minus) the probability current. As a function of the 
initial condition Uj, Qt{u) = V{uf; tf\ui = u;ti = t) satis- 
fies the backward Fokker-Planck equation 

-dtQt{u) ^ udlQtiu) - {ii- wt)duQt{u) 



dluQt{u) +du{wt 



u)Qt{u) + Qt{u). 
(122) 



The propagator also satisfies the initial condition 

Viuf, ti|ui; ti) = V{u[; tt|ui; tf) = S{u[ - Ui). (123) 

In order to obtain the solution of the forward or backward 
FPEs, one needs to complement them with a boundary 
condition at w = 0. We consider two cases: 

(i) Propagator with reflecting boundary Prefl- This is 
the case we studied so far in this article, and in [JS]. It is 
defined by a vanishing probability current: Jt(u = 0+) = 
0, or 



lim duuPt{u) + (u- wt)Pt{u) = 0. 



(124) 



Typically, this is satisfied by a power-law-like behaviour 
P{u) ~ w^i+"'* for u ^ (see the examples in [20l [24l 

m)- 

(ii) Propagator with absorbing boundary Vahs- This 
is the relevant case for sub- avalanches. The problem is 
characterized by a vanishing propagator, when starting 
from ii — 0~^, i.e. Qt{u = 0+) = 0. This implies (except 
for pathological cases) that the "current" for the back- 
wards FPE vanishes. 



lim duuQt{u) -t- (w* - 2 - u)Qt{u) = 0. 

ii— fO 



(125) 
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On the other hand, for the forward Fokker-Planck equa- 
tion, Pt{u = 0+) = Prediui = 0+,if = t\ui;ti) will typ- 
ically be a non-vanishing, non-trivial function of time, 
and the current Jt will not vanish (as expected from 
physical intuition, since trajectories touching ii = are 
"absorbed" ) . This is why treating the absorbing bound- 
ary using the forward Fokker-Planck equation is inconve- 
nient; instead, the backward equation is natural hercFM 
Let us now define Laplace transforms with respect to 
the final velocity for P and with respect to the initial 
velocity for Q: 



Pt{X\uuti) 



duie^'^'V{ui;tf^t\ui;ti), (126) 



Qt{X\uf, tt) := / dui e^"'7'(wf ; Ulin; U = t). (127) 
Jo 



Laplace-transforming the forward Fokker-Planck equa- 
tion (120), with the reflecting boundary condition (124), 



we obtain a first-order PDE for Pf(A) 

dtPr^'W = (A^ - A) ^A^'^^A) + u;tAPr^A). (128) 

Note that while boundary terms would arise in gen- 
eral, the reflecting boundary condition ( 124 ) ensures that 
they vanish. Similarly, Laplace-transforming the back- 
ward Fokker-Planck equation (122) with the absorbing 



boundary condition (1251, wc obtain a first-order PDE 
for Qt{\): 



-<9,Qf ^(A) = (A^ 



X)dxQt'{\) 

+ [{2-wt)X + l]Q 



(A). 



(129) 



Again, vanishing of boundary terms for the Laplace 
transformation is ensured by the absorbing boundary 
condition (125). On the other hand, the Laplace- 



transformed equations for P^ and Q\ ^^^ more com- 
plicated: There, the boundary terms at m = do not 
vanish and are undetermined functions of time. Thus, in 
the following, we will always use the propagator in terms 
of the final condition P or its Laplace-transform P when 
discussing a reflecting boundary, and the propagator in 
terms of the initial condition Q or its Laplace-transform 
Q when discussing an absorbing boundary. 

Eqs. ( 128 ) and ( 129 ) can be solved using the 



Now, 



method of characteristics. In the forward (reflecting 



^^ Similarly, for the reflecting boundary, Qt(u = 0+) = 
'Pi^i\{ui,tf\ui = 0+;ti = i) will typically be a non-vanishing, 
non-trivial function of time. So, for the reflecting boundary, the 
backward equation is inconvenient and the forward equation is 
natural. For the absorbing boundary it is the other way around. 
This peculiar behaviour is due to the nature of the ABBM noise 
term, which vanishes for it = 0. For a standard Brownian mo- 
tion, the absorbing boundary can be treated equally well using 
the forward or the backward Fokker-Planck equation: There we 
have "PabsC^f = 0) = Pabsl"! = 0) = 0. 



boundary) case this method was shown to provide a gen- 
eral connection between the Fokker-Planck approach to 
the ABBM model and the dynamical path integral (in- 
stanton equation) approach [T51 US] . In the following we 
take a flrst step towards generalizing this to the case of 
an absorbing boundary. The solution of ( 129 ) is 



Q,,(u(iO)=e^''[(^--)-+il'^*Q,,(^(if)), 
where it satisfles the backward instanton equation 
- dtiit + ut + ill = 0, 



(130) 



(131) 

with the boundary condition ii{ti) = AFj 

The initial condition (123) for the propagator gives 
Qtt{Xf) = e^f^^ . Inserting this into (130), we obtain 



ditie^'"7'abs("f;tf|"i;ii) = Qti(A) 



exp 



[(2 - Ws)Us + 1] ds + UfMtf 



(132) 



It is useful to recall that the fact that this solves the 
problem with an absorbing boundary is an indirect conse- 
quence of our chain of arguments: It stems from the fact 
that we use the backward instanton equation ( 131[ ) which 
encodes the solution (via the method of characteristics) 
of the Laplace-transformed backward FPE, which itself 
contains no boundary term precisely in the case of an ab- 
sorbing boundary. For the case of a reflecting boundary, 
the analogous formula is 



duf e^"f7'rcfi(uf;if|ui;ti) = Af(A) 



exp 



Ws'u(s) ds + iiiu{ti) 



(133) 



As discussed in Ref. [4<F , this is equivalent to solving ( 19 ) 
as above using the MSR field theory and the instanton 
equations. For details, see [IHj section V C, and, in the 



present article, section VIII D where we discuss this in 
general for the ABBM model with retardation. 



1. Propagator at constant driving 



Let US now apply ( 132 ) in order to determine the prop- 



agator of the standard ABBM model with an absorbing 
boundary, at a constant driving velocity < t; < 1. 
The solution of Eq. (tml) is 



u{t) 



xe{t-ti) 

s*i-*(A + l)-A' 



(134) 



^^ Note that this equation is causal — solved by increasing the time 
— in contrast to the forward instanton. 
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Inserting this into ( 132 ), we obtain the Laplace transform 
(w.r.t. the initial condition) of the propagator with an 
absorbing boundary, at a constant driving velocity v, 



duie^"'7'abs(Mf;if|Mi;ii) = Qti(A) 



cxp 



cxp 



[{2-v)u{t) + l]dt + u{t{)u{ 



ti 



ti 



{2~v) / u{t)At + [ti - ti) + u{ti)ui 

Ae*f"*'Uf 
^- \- ft ^ f- 

1 - A (e*'-*^ - 1) ' '. 



exp 

X [1 - A (e*'-*' - 1)] 



(135) 



Inverting the Laplace transform from A to iii yields 
the propagator of the ABBM model with an absorbing 
boundary at u = 0, 



7^abs(wf;tf|ui;0) = 



(136) 



exp -if 



-Ufe* 



1 



jUf/Uj) ^ 

2 sinh I 



-Il-v 



V sinh I 



Here we set ti — for simplicity since the result depends 
only on tf — ti. 

Our result is identical to Eq. (37) in the recent cal- 
culation [SB], there obtained using completely different 
methods (decomposition in eigenfunctions of the Fokker- 
Planck operator) . The advantage of our approach is that 
it makes the connection to field theory clearer, and that 
it is easily generalizable to situations with a non-constant 
driving velocity (where the eigenfunction method is not 
applicable) . 



It is straightforward to check explicitly that ( 136 ) sat- 
isfies the backward FPE (122). Near Ui = we have 
7'abs(wf; tf|ui; 0) ^ ul~". For < w < 1, it satisfies the 
absorbing boundary condition T^abs (wf ; if | Ui = 0;ii) = 
and ( 125 ). On the other hand, near iif — 0, T^abs is a non- 



vanishing constant, so the "forward" current (121) does 
not vanish. For velocities u > 1, our assumption on the 
absorbing boundary condition, 7^abs(uf; tf |ui = 0;ii) = 
and (125), seems to break down. In fact, the absorbing 



boundary becomes equal to the reflecting one since ii = 
is unreachable!^ 



2. Sub-avalanche durations m the standard ABBM model 

Another simple example is the "survival probability" , 
i.e. the probability never to touch the absorbing bound- 
ary u = until time tf , starting at some initial iii at time 



^^ In fact, for i^ > 1, one can still formally define avalanche-size 
distributions conditioned to ii being close to instead of precisely 
as for <.v < 1, see the appendix in I45| . 



ti-. 



^surv(tf|Wi;ti) 



AuiVs.hs{ui]ti\ui\ti) 



Its leading behaviour as ui — )■ corresponds to the prob- 
ability of an avalanche duration T > tf — ti. 



It can be obtained using ( 130 ) as in the previous sec- 



tion, but with the initial condition 



QtAh) - 



dii t 



1 



(137) 



for A/ < 0. 

As above, inserting the solution (134) for u(f), we ob- 
tain for the survival probability at time if 



fdUie^'''Psurv{tf\Ui;ti) 



exp 



[(2-u;^)u,, + l]ds 



1 



[1 + A - Ae*f~*'] 

[l + A-Ae*^-*f 
A 



U(if) 

v-2 (A + 1)-Ae*^-*- 
A 



(138) 



This is, of course, identical to the integral of ( 135 ) over 



iif. Inverting the Laplace transform, we obtain the sur- 
vival probability at time if > ij when starting from -u; at 
ti . It is a function of if — ij only. For w < 1 it reads 



-* survl^^f ''il^ij 



1 



r(i-^) 



(139) 



It increases from to 1 as u increases from to oo, while 
for u > 1 it is equal to unity for all u > 0, i.e. there 
are no zeroes of the velocity. Note that this result was 
obtained independently in |59j and |58j by completely 
different methods. It can also be obtained by integrating 
the propagator (136) over iif from to oo. 



The case w < 1 is considered from now on. Taking 
a derivative w.r.t. the final time gives the probability 
density of first-passage times Tq for ii to become zero, 
given an initial velocity Wi, 



-PfirstCTolUi) 



r(l-„)(eTo_i)2 



e^O-l 



+To 



(140) 



Since an avalanche always starts at u = 0^, we can ex- 
tract from the leading term in small iii a density of du- 
rations Tq, 

,^ ■. T , s / -r nu-2 sinfTTt;) 

«Pduratio„ To = e^« 1 - v) e^« - 1 -^-^ . 141 

In the limit w — )- 0, density means units of 1/u = 1/w. 
Since /5duration(T'o) diverges like 1/Tq~ , hence is not nor- 
malizable for v < 1, we have chosen to normalized it as 
{vTq) — 1. Note that p has a finite limit at w = 0"*" which 
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agrees with the avalanche-duration density obtained in 

The probabihty density for an avalanche to continue a 
time To beyond an arbitrary chosen time ii is obtained by 
integrating over the distribution of uq in the stationary 

state Pstat(Mo) = ""o ^^"tcSJ"' 



-Pbeyond(To) = / dUj F^tat ('ii)-Pfirst (To |ui) 

Jo 



(142) 



Note that this is a bona-fide probability distribution (nor- 
malized to unity). 

Finally, we obtain again the density of avalanche du- 
rations by taking a derivative, 



vp{T) 



-d, 



(e^ 



I T=To 

T 
f 



-'beyond \-^) 

sin(7rf ) 



(143) 



This density can be interpreted as a probability vp{T) = 
P{T) in the following sense: Take a random time t — 0. 
The velocity uq is positive with probability 1. Thus, there 
is a first-passage time —ti at zero velocity on the left of 
t = 0, and a first-passage time +^2 at zero velocity on the 
right of t = 0. The duration T = t2 - (-ti) = ^2 + ti is 
the distance between the two. It is thus normalized not 
by /g°° P{T)dT = 1, but rather by 



dii 



At2P{T = t2+h) = l. 



In terms of units, P{T) is not a probability density 
(which would give a dimensionless number when multi- 
plied by a time interval dT), but a "double" -probability 
density which gives a dimensionless number when mul- 
tiplied by two time intervals dii, di2- In other words 
J TP{T)dT = 1 since the probability that a randomly 
chosen time belongs to an avalanche is proportional to 
its duration T. 

Note that for w — >■ Eq. (143) reduces to the result 
known from [211 113 HI] ■ The velocity-dependent power 
law for small T, P{T) ex T~^+^, was already predicted 
in [23j . A similar result for the distribution of avalanche 
sizes at finite velocity in the standard ABBM model is 
discussed in appendix IE) 



B. Fokker-Planck equations and propagator 
including retardation 

Now let us go back to the more general ABBM model 
with retardation. 



Fokker-Planck equatiorp^ generalizing (144 1 for the joint 
probability distribution P{u, h), 

dtPt{u,h) 

= d^uPt{u, h) + du [u — Wt + aii — ah) Pt{u, h) 
+ T-^dh{h-u)Pt{u,h). (144) 

As discussed in [1^, this forward Fokker-Planck equa- 
tion provides an alternative derivation for the generat- 
ing function (26). The instanton equations (24), (25) 



are equivalent to the equations for the characteristics of 
the linear PDE ( |144[ ), see [49, section V C for details. 
The transformation between the "real space" u, h and 
the "Laplace space" m, h of the characteristics, i.e. in- 
stantons is a very useful tool whenever boundary terms 
are absent. This is the case for a zero probability current 
at u = 0, i.e. for a reflecting boundary condition. 

To study the case of an absorbing boundary, as noted 
above for the pure ABBM model, it is useful to consider 
the flow of the probability density as a function of the 
initial condition u, h at time ii = i, which satisfies the 
backward Fokker-Planck equatiorP^l [gp] 

-dtQt{u-,h) 

= ud\Qt{u, h) — [ii ~ Wt + ail ~ ah) duQt{u: h) 

-T-^ih-u)dhQt{u,h) (145) 

= d^uQt{u, h) — du {ii + 2 — Wt + ail — ah) Qt{u, h) 

- T-^dh (h - it) Qt{u, h) + {l + a + T-^) Qt{u, h) 



Both ( 144 ) and ( |145[ ) are linear in the probability density 
P 01 Q. Hence, they are completely characterized by 
the Green function or propagator 7'(-Uf , /if ; ijuj, /i;, 0), the 
probabili ty to go from iii,hi at tj to iii,hi at ti > ti. It 
satisfies ( 144) as a function oi t = tf^ii = ii,[,h — hi and 
(1451 as a function oi t = tf,ii — iii,h — hi] it has the 
initial condition V{iif,hf;ti\iii,hi,ti) — (5(uf — Ui)S{hf — 
hi). 



To derive the forward Fokker-Planck equation | |144| |, we set 

Ptiu,h) = (5{u-ut)5ih-ht)) . 

Then apply Ito calculus to 

Pt+dt{u, h) = {5{u - [iit + dut])5{h - [ht + dht])) , 

with 

diit = ViitdBt + [{wt - iit) - aiit + aht]dt , 

rdht = [iit — ht]dt 
(dBtdBt,) =25{t-t')dt. 



For the backward equation ([145|, apply Ito calculus to 

Qt+dt{ut+dt,ht+dt) = Qt{ut,ht) . 



The equations of motion ( 23 ) are equivalent to a 
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Monotonicity, domains of definition, and 
boundaries 



It is important to note that ii and h satisfy together a 
monotonicity property: If li; > and both u{t = 0) > 
and h{t = 0) > 0, then they remain so at aU times. Al- 
though the quadrant m > 0, /i > is the more physical 
one, as we see below it is convenient to solve the FP equa- 
tions in the half-space fi > and then at the end restrict 
to the quadrant /i > 0. The reason for this is simple: No 
matter whether we impose an absorbing or a reflecting 
boundary at ?i = 0, one has a finite probability of reach- 
ing Mf > 0, /if > starting from Uj > 0, /li < 0. So, as 
a function of the initial value /i;, 7^(uf, /if;ii|ui, /ij, tj) is 
smooth around hi = 0, and has a finite value there. Its 
natural boundary is at hi — >■ — oo, where V ^ 0. Thus, 
when considering the backward equation, we will work 
on the half-space Mi > 0, ft-i € M in the following, in order 
to avoid undetermined boundary terms. 

Due to the aforementioned monotonicity property, we 
know that the restriction to /ij > of the half-space prop- 
agator obtained in this way, will actually be equal to the 
propagator restricted to the physical quadrant from the 
beginning. 

By taking linear functionals of the propagator, one ob- 
tains other observables. We give a few examples: 

(i) Starting not from a fixed point, but from a distri- 
bution of initial values Pi{ui, hi), the probability density 
of final values 



Pt{u,h) 



dui 



dhi Pi{ui, hi)'P{u, h; t\ui, hi; ti) 



(146) 
still satisfies the forward FPE ( 144 1 . The initial condition 
is, naturally, Pt- = Pi- 

(ii) One may be interested not in the probability den- 
sity of the final point at U{, h{ (given by the propagator), 
but the actual probability to land in a domain D, start- 
ing from an arbitrary initial condition u, h. This is given 
by 

Qtiu,h)= (147) 

) /"OO 

d-Uf / dhiQi{uf,hf)'P(uf,h{;t\u,h;ti), 



where Qi{uf,hi) = Id is 1 inside D and outside. 



Qt{u,h) satisfies the backward FPE (1451, with the ini- 
tial condition Qt^ ^ Qi = '^d- This can be used to de- 
termine the distribution of avalanche durations, starting 
from an initial value Uj > 0, /li. Choosing D to be the set 
{iii > 0}, one obtains the probability to have any posi- 
tive domain- wall velocity at if, i.e. the probability not to 
have touched the boundary u = between ti and U 



As mentioned above, our instanton solution (26 1 is. 



equivalent to the forward propagator V{uf, hf;t{\ui,hi,ti) 
from a given initial condition iii — 0,hi — to some fi- 
nal point ill > 0, /if > 0, with a reflecting boundary 
at u = (the boundary at /i = is unreachable when 



propagating forward). Imposing an absorbing bound- 
ary at u = 0, as required for analyzing sub-avalanche 
durations, is less trivial. In contrast to the case of 
e.g. a standard Brownian motion, the probability cur- 
rent at the final point iif ~ vanishes as soon as one 
sets Vahsiui — 0,h{;t\ui,hi;ti) = 0. The correct prop- 
agator with an absorbing boundary should thus have 
'Pabs(uf = 0,h{;t\ui,hi;ti) > 0, an undetermined func- 
tion of time, as confirmed by the explicit calculation in 
section VIII A 1 (see also the discussion in Appendix D 



of ,18J). Hence obtaining it from the forward FPE (144) 
is not easy. However, in terms of the initial condition, as 
motivated for the pure ABBM model above, we expect 
Vahsj ut, hi] t\ui == 0, /li, ti) = 0. Then, the backward FPE 



( 145 1 is easy to analyze using Laplace transforms, since 
the boundary term Q{ui = 0, hi) vanishes. 



D. Characteristics and instantons 

To solve the backward FPE, we define the Laplace 
transform Q via 



U\^^)■■= du d/ie^"+^'^Qt(u,/i). (148) 

Jo J-oo 



Equation ( 145| then gives 
~dtQtiX,tJ-) = [A" + (1 + a)X - r-i/^] DxQtiX, m) 



[-aX 



V] df,QtiX,n) 



+ [i2~wt)X+{l + a + T-i)] 4(A, n). 

(149) 

Note that to obtain this equation (with vanishing bound- 
ary terms) we used Qt{u = 0,/i) = 0, and the fact that 
the noise vanishes at m = (which is a special property 
of the ABBM model). [3 

We now define the characteristics u, h via the backward 
instanton equations 

-dtu{t) + (1 + a)u{t) + u{t)^ - h{t) = (150) 

~Tdth{t) + h{t) - au{t) = (151) 

~dtq{t) = [(2 - wt)u{t) + (1 + a + r-i)] q{t) (152) 



1^ Eq. |T49l is valid for the LT | |148| defined on the half-space u > 
0, /i £ M. If one wishes to define the LT on the physical quadrant 
u > 0,h > 0, one cannot avoid a boundary term from h = 0, since 
Qt{u,h = 0+) does not vanish. Since there is no noise term in 
the evolution equation for h, the solution does not have to be 
continuous, and there is no contradiction to Qt{u,h = 0) = 0. 
One possibility to eliminate the boundar y ter m is to add an extra 
diffusion term ehd^Qt on the r.h.s. of l |145[ l, which leads to an 

additional te rm e^ ^df^Qt on the r.h.s. of l |149| l, and eh'^ on the 
l.h.s. of Eq. ( |151[ l. Then Qt vanishes at h = 0, is continuous, 
and one can show that as e — > it converges pointwise for any 
h > to the (discontinuous) restriction to the quadrant of the 
presently considered solution in full space. Wc will however not 
need to further explore this method here, but can work with h 
on the full real line. 
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The boundary conditions are u{ti) — A, h{ti) — fi/r, 
q{ti) = Qti(A,/x). Note that the first two equations are 
identical to the standard instanton equations (24 1 and 



(25 ). The equation for Q{\, /i) along a characteristic sim- 
plifies to 

- -Qtiu{t),Th{t)) 

= [i2-wt)uit) + il + a + T-^)]Qt{u{t),Thit)). (153) 
Its solution is 

Qt,{u{ti),Th{ti)) 

= exp i / [(2 - Ws)f4s) + (1 + a + r"^)] d.s I x 

xQt,{u{ti),Th{tf)). (154) 

The initial condition Qt^ depends on the observable we 
want to compute. For concreteness, let us think about 
the propagator T^abs- There we have 



Qtf {u, h) —Siii — U{)S{h — hi) 



Qtf (u{ti),Th{ti) 



^uiu{ti}+Thth{ti) 



(155) 



Then, the propagator Pabs is given by inverting the 
Laplace transform 



dui 



.(A,m) 



exp 



dhi e^"'+^'*'7'abs(uf , hf, tf\uu h; U) = 



[{2-^Ws)u{s) + {l + a + T-^)]ds (156) 



+Uiu{tf) + Th{h{tf) > 



To summarize, one can say the following: The usual 
instanton equations (24) and (p5|) encode the solution 



of our model with a reflecting boundary at w = 0. 
The backward instanton equations ( 150 ) and ( 151 ) en- 



code the solution with an absorbing boundary at w = 
0, assuming the absorbing boundary indeed satisfies 
'Pa.hs{ui,h[;t\ui = 0+,hi,ti) = 0. 

It would be interesting to extend this approach to a 
general retardation kernel, beyond the simple exponen- 
tial. In that case no local-in-time Fokker-Planck ap- 
proach seems available. It is tempting to conjecture 
that the correct generalization of the backward instan- 
ton equation ( 150 1 is 



-dtut + (1 + a)ut +u'f + a I ds f'{t - s)us = -At, 

(157) 



corresponding to ( 18 ) after mapping u — >■ —u. Note that 
( 157[) r educes to ( 150 ) and ( 151 ) for the exponential ker- 
iieTpl), by setting 



^{t-s)/T 



iisds, 



(158) 



as would be required. 



E. Sub-avalanche durations with retardation 

In order to compute sub-avalanche durations, we need 
to eliminate the h variable from (156). This is done by 
integrating over hf and fixing a value of /liFj In appendix 
[F| we motivate the following generalization of ( 156 1 for 
this case. 



oo poo 

dhf / dui e^"' Pabs (wf , /if ; if I Mi , /li ; ii ) = 

-oo Jo 

cxpj / [(2- Ws)us + (1 + a + T"^)] ds 
+UfUtt — rhiht- > X 



8- 

'-'hitt)\h(tt)=0 



Mu) 



(159) 



Now u, h are solutions of the backward instanton equa- 
tions (150), (|151[) with the boundary conditions 



h{ti) = 0, 



u{ti) = A. 



(160) 



The effect of going from a fixed hi to a fixed hi is thus 
a change in boundary conditions for h, and the Jacobian 



factor in (159). To understand its importance, let us see 
how (159) reduces to the pure ABBM solution (132 1 in 



the case of a = 0. This is not trivial, for example the 
exponential factors in (156), (159) and (132) are quite 
different. 



1. Recovering the pure ABBM model 



For a = 0, we can solve Eq. (151) explicitly. 



hi = /iie(*'-*')/^ 



(161) 



Thus, the Jacobian factor in (159) is 



d, 



h{U)\hitt)-- 



__,hiU 



.) - f,-itl-ti)/r 



Inserting this into ( 159 ) for a = gives 



dh, 



f / duie^"''Pi,hs{uf,hi;t[\uuhi-ti) = 

-oo Jo 

exp < [(2 - Ws)us + 1] ds + ?ifUtj - rhiht, > , 

(162) 



where now h, u are solutions of (150) 



and with the boundary conditions ( 160 



TSOl), jTSl] ) with a 
s ( |160p . Note that 



= 

the 



^^ Or convoluting later with a normalized distribution of h^. In any 
case, this is not the same as setting /i = in l |156[ l. The latter 
would be an integral over all h;, which in general diverges. For 
example, in the pure ABBM model the sub-avalanche duration 
is independent of /i;. 
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Jacobian factor was necessary to cancel the term (if 



ii)/r of the exponential. The solution (161) for h with 



the boundary condition hf — from ( 160 1 implies that 



ht — for all t. Hence, ( 162 ) reduces to the pure ABBM 
solution (132), and the equation (150) for u reduces to 



the pure ABBM instanton equation ( 131 ) 



As expected, we obtain that in the limit a — 



dhf'Pi,hsiuf,hf]t{\ui,hi;ti) = Vs,bsiut;tf\ui;ti) 



is independent of h^ and reduces to the propagator of the 
standard ABBM model. 

Now let us return to the more interesting case of the 
ABBM model with retardation. In the following section 



we will apply (159) to determine the correction to the 



duration distribution of the first sub-avalanche, for small 
a (weak relaxation) and r = 1. This is similar to the 
perturbation theory in a performed in appendix [Cj We 
also attempted to analyze the backward instanton solu- 
tion in the physically interesting limits of fast and slow 
relaxation, as done in sections IV] and VI for the forward 
(reflecting boundary) solution. However, we encountered 



technical difficulties and leave this for future work. 



2. Weak-relaxation limit 



Let us consider the solution of ( 150 ), ( 151 ) with bound- 
ary condition 



Hti) 



M, 



u{ti) = A. 



(163) 



We will be interested in the limit of small fi (required for 
computing the Jacobian factor in (1591; the rest can be 



computed at /Lt = 0) . To this end we expand the instanton 
to order a, ^ as 



fik 



(10) 



ah. 



(01) 



+ afih['^^ +0{a^,fj,^). 



(164) 



The correction of order a to Psurv, i-C- the cumulative 
distribution function for the avalanche durations, is ob- 



tained from (1591 as 



-* surv '^'^i "-i j ■- 



exp 



= exp 



duie^"'Psurv(if|"i,/ii;^i) 



/•OO /"OO 

duie^"' / diii 

Jo J-oo 



[(2 - Ws)us + {l + a + T ^)]ds- Th-M, 



1 



u{U) 



dhiVahs{uf,hf\tf\ui,hi]ti) 



(2 - w,) («(°") + auT^) + (1 + a + r'^)] ds - arhhf^^ 



P^=':{\)ll + a 



(2 - ws)uf^'^ + 1 



ds 



,^,(01) ^(11) 



^(00) f (10) 



1 



exp 



ii("o)(tf)-^au(oi)(if 



(01), 



-aTh). hi 



K 



(10) 



;;(11) 
ahl ' 



(165) 



The pure ABBM result P°=°(A) is given by ( pSf . By 
time translation invariance, this expression only depends 
on T := if — t\. In appendix [D] we perform the pertur- 
bative calculation and determine explicitly the functions 
appearing in the expansion (164 1, for the case r = 1. In 



figure [15^ we show the form of the resulting backward in- 
stantons, obtained from the pcrturbative expansion and 
a numerical solution. The case of general r is more com- 
plicated and left for future research. 

The survival probability Psurv is the probability of hav- 
ing a first-passage time to m = of T > if — ij. The (nor- 
malized) probability distribution of first-passage times at 
w = 0, starting from mj is thus 

Pfirst(T'|Ui, /li) = 9tf |j^.^j.^j,Psurv(if|Ui, /ii;ii). (166) 



By time-translation invariance, this is independent of i\. 
For the pure ABBM model, we computed Pfirst in ( |140| ). 
The (unnormalized) density of avalanche durations jO(r) 
at a fixed value of h\ is the leading order of -Pfirst (wi, h\) 
as Mi — >■ 0, in our case 



p{T) :^ lim w°'''+" ^ Pii^stiT\ui,hi). 

iii— ^0 



We can then express p{T) as the leading order of Pg, 

as A — > — oo: 
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p(T')cir'= lim P(A,0)(-A) 

A— >- — CX3 



-2+t»+a/ii _ 



(167) 



-ahiT 



^ ^ ^ e'^[~2{v^l)U,{l-e'^)-{T' + l)v + T' + T-l]+{2-v)T + v + l ^ ^^^^^ 



(^eT - i)i-v-ah, I ' 2{eT -1) 

This further simpUfics for the case of a vanishing driving velocity, w — >■ 0. There we obtain 



p{T')dr 



fpT _ l\l-ahi 



1 + a- 



eT (t2 + 2Li2 (1 - e^ 



T- 1 +2T+1 



2(e2 



1 



Oia)' 



On the other hand, expanding (167) for small T we get 



0.010 



ahi + V ~ 1 



(168) 



+ hah, + v)il + a) (v - 1 - ^ ] T 
2 V 1 + a. 



©(T)^ 



We observe that the power-law behavior near T = 0, 
p{T) ~ y-z+w^ jg j^Q^ modified for the first sub-avalanche 
{hi = 0), but is modified for later ones to 



p{T) ^ T 



-2+v+ahi 



(169) 



to leading order in a. This is natural, since for small 
avalanches h remains essentially unchanged, and replaces 
in ( 19 ) w — > v + ahi. In order to obtain the sub-avalanche- 
duration distribution for stationary driving, one would 
need to average over all hi taken from the stationary dis- 
tribution P{hi\ui = 0). Presumably, this would replace 
the correction ahi to the exponent in ( 169 ) by a velocity- 
dependent correction. We leave the details for further 
research. 



In figures I9] and 10 we compare the result (167) to nu- 
merical simulations of the original model. One observes 
good agreement for small a but larger deviations starting 
around a w 0.5. In figure 11 we verify numerically the 



result ( 169 1 for the sub-avalanche-duration exponent as 
a function of /i;. The agreement is very good, even for 
a = l. 



3. Numerical results 

Since the analytical results we obtained above for sub- 
avalanche sizes and durations are rather limited, we also 
give a few qualitative numerical results. In this section we 
neglect the difference between densities and probability 
distributions, and use P instead of p everywhere. 

1. Sub-avalanches at fixed initial Ui^ hi — Q and v — 0. 
In (168 1, we derived that the ABB M power-law 
P{TY^ r-2 for smafi T (see plsj )) remains un- 
changed in this case, at least for t — 1 and small 
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Figure 9. (Color online). Duration distribution of an 
avalanche starting at u = 0, ^ = and t = l,v = 0.6. 
Crosses, diamonds and circles: Results from numerical in- 
tegration of ( |19[ ) for a — 0, a = 0.3, a = 0.5 with time step 
di = 10"''. Lines: Expansion in a, (1671. Yellow (thick) line: 
Pure ABBM model, a = 0, Red (dashed) line: a = 0.3, Black 
(dotted) line: a — 0.5. 



a. In figure 12 we consider the sub- avalanche 
sizes and durations for general values of a,T. We 
see that the mean-field zero-velocity power-laws 
P{S) ~ S-^/'^,P{T) ~ T-2 are clearly visible even 
when varying a, t. For large a and small r, there is 
an interesting crossover in the shape of the distri- 
butions, showing a similar power law but with dif- 
ferent amplitudes, depending on a,T. The case of 
large t is equivalent to a modification of the mass, 
i.e. of the large- avalanche cutoff, as discussed in 
section III Al 

2. Sub- avalanches at fixed iii, and hi taken from the 
stationary distribution at the driving velocity. We 
observe in figure [13] that this leads to a modifica- 
tion of the ABBM power-law exponent. However, 
the variation in the exponent due to retardation 
becomes smaller as the driving velocity decreases. 
This is expected from (169), since the typical value 



of hi (and hence the correction to the exponent) 
decreases as z; — )■ 0. 
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Figure 10. (Color online). Correction to the duration distri- 
bution of an avalanche starting at u = 0, ^ = due to re- 
tardation, [P°(r) - P''=°(T)] / [aP^^^iT)]. r =l,v = 0.6. 
(Yellow) circles, (red) squares: Results from numerical inte- 
gration of ( 19 1 for a — 0.1, a = 0.3 with time step dt = 10~*. 
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Figure 11. (Color online). Duration distribution of the first 
sub-avalanche, as a function of /li. r = 1, a = 1, u = 0. 
(Yellow) crosses, (Red) diamonds, (Blue) circles and (grey) 
squares correspond to numerical simulations of ( |19[ l for hi = 
0, hi = 0.1, hi = 0.2, hi = 0.3 with time step dt = 10"^. Solid, 
dashed, dot-dashed and dotted lines correspond to power laws 



From these numerical results we conjecture that sub- 
avalanche durations in the ABBM model vifith retarda- 
tion, at constant driving, satisfy 

P(T) r^ y-2 + y(l+CT) pfg\ ^ g-3/2+v/2{l+cs) 

Here ct, cs depend on a,T, and vanish as a — >■ and/or 
T — > 0. In other words, we conjecture that the zero- 
velocity exponents are unchanged, and only the driving- 
velocity-dependent part is modified. This is also con- 
sistent with our analytical results in sections |V] and VI 
for the velocity distribution P{u). Verifying these con- 
jectures in more detail, numerically or analytically, is an 
interesting task for the future. 



IX. CONCLUSION AND OUTLOOK 

In this article we have analyzed in detail a general 
ABBM model with retardation. We showed that it sat- 
isfies the Middleton property (monotonicity of the dy- 
namics). Using this, and the Brownian nature of the 
ABBM disorder, we have been able to reduce the calcu- 
lation of the expectation value of a general observable 
in presence of monotonous driving to the problem of 
solving retarded non-linear instanton equations. These 
equations can be implemented numerically for arbitrary 
retardation kernels, and are much simpler than the orig- 
inal model with quenched disorder. To obtain analytical 
results, we focused on a model with exponential relax- 
ation, which reduces to two coupled "instanton" equa- 
tions, local in time. We have derived explicit forms for a 
number of observables at stationary and non-stationary 
driving. We mostly focused on the two limits of fast 
relaxation (describing eddy current effects in magnetic 
Barkhausen noise) and slow relaxation (of interest for 
earthquake models). 

In the limit of slow relaxation, the main physics is 
the splitting of a single avalanche of size S of the stan- 
dard ABBM model, into a "cluster of aftershocks" , of the 
same total size S — "^^ Sa, but of much longer duration 
(strictly infinite for an exponential retardation kernel). 
This splitting is sharp in the limit of quasi-static driv- 
ing and strong time-scale separation. Although we have 
been able to quantify some of these ideas, a more detailed 
analysis of the statistics of these sub- avalanches remains 
an important challenge for future work. In particular, if 
these ideas are to be extended to realistic earthquake dy- 
namics, one needs to investigate power-law retardation 
kernels motivated by the Omori law of decay of activ- 
ity. Still, it is of great interest to have found a tractable 
model with aftershocks, given that the standard ABBM 
model, which also models interfaces within mean-field 
theory (at the upper critical dimension) yields indepen- 
dent avalanches following a Levy process [T7| [T51 [53] |_] 

In the context of magnetic systems, the retardation 
describes the influence of eddy currents on the statistics 
of Barkhausen noise pulses. Experimentally, the time 
scale of eddy current relaxation is much shorter than that 
of the domain-wall motion, motivating our study of the 
fast-relaxation limit. In this limit, the effects of retarda- 
tion are perturbative and vanish when the eddy-current 
relaxation-time tends to zero. We have obtained the lead- 
ing corrections to the stationary velocity distribution, as 
well as the decay following a kick in the driving velocity. 

In both the slow-relaxation and the fast-relaxation 
limit, the influence of retardation turned out to be most 
pronounced in non-stationary observables. For example. 



^^ While this work was completed we learned of an independent 
work by Jagla, Rosso |61| . These authors consider several vari- 
ants of earthquake models, and one of them is similar to the 
retarded model considered here. 
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Figure 12. (Color online) First sub-avalanche sizes (right) and durations (left) from numerical simulation of the ABBM model 
with retardation, starting from iii = 0, fti = 0, at u = 0. Grey (dashed) lines: power laws P{S) ~ S~^'^, P{T) ~ T~^. Yellow 
(thick) line: pure ABBM model a = 0. Red (dashed) line: a — 60, r — 0.05. Blue (thin) line: a = 5,r = 10. 
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Figure 13. (Color online) Sub-avalanche sizes from numerical simulation of the ABBM model with retardation, in a steady 
state with v = 0.3 (left) and v = 0.02 (right). Yellow (dashed) lines: power laws P{S) ~ S-ta-'')/^^ left: S'"^■^^ right: S"^*^ 
Red (dotted) lines: power laws, left: S~^'^ , right: S~^'^^. Crosses: pure ABBM model a = 0. Diamonds: a = 1, Circles: a — 2. 
In all cases, r = 1. 



we computed explicitly the tail of the avalanche shape at 
fixed size. While it is Gaussian in the standard ABBM 
model, in the ABBM model with retardation it decays ex- 
ponentially. Furthermore, we showed that formally the 
avalanche activity following a kick in the driving veloc- 
ity never stops (even if in the fast-relaxation limit, all 
significant avalanche activity still ceases), also in strong 
contrast to the pure ABBM model. 

One important direction for future theoretical inves- 
tigations is the role of the internal dimension of the in- 
terface. Here we reduce the description of the magnetic 
domain wall, which is a two-dimensional elastic interface 
in a three-dimensional medium, to a single degree of free- 
dom (its center of mass) . This mean- field description has 
been argued [21] and recently shown in detail [iSl US] to 
be accurate for the center of mass of the interface above 
a critical internal dimension dc- For certain soft mag- 
nets, due to long-range elastic forces, one has d^ = 2 
|24j . Hence, some realistic domain walls have just the 



critical internal dimension. In that situation the inter- 
nal degrees of freedom of the interface contribute only 
logarithmic corrections to the mean-field behavior as de- 
scribed in detail in [18j . Other ferromagnets are known to 
be in a non-mean-field universality class |31j . There, cor- 
rectly describing the details of Barkhausen noise requires 
combining the eddy current modifications discussed in 
this article with a treatment of the spatial degrees of 
freedom, for example using the functional renormaliza- 
tion group as in [TSl [551 112] ■ 

Another important avenue for further work is a de- 
tailed comparison of the analytical results discussed here 
with experiments on real magnetic domain walls. Consid- 
ering the above results, non-stationary observables seem 
most promising. For example, the avalanche shape dis- 



cussed in section VII B shows a clear qualitative differ- 
ence to the result of the standard ABBM model. As 
recognized in [121 HI] this is an inertial effect due to an 
effective negative domain-wall mass. However, inertia 
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can be modelled in different ways - for example, the re- 
tarded ABBM model considered here and in [22], the 
ABBM model with second-order dynamics considered in 
[15] , and the mean-field model with stress overshoots con- 
sidered in [10]. We believe the avalanche shape allows 
identifying not just the existence and the sign of iner- 
tial effects, but their precise form. To this end, more 
precise analytical and experimental results (which go be- 
yond a single "skewness" quantity, and consider the en- 
tire shape) are necessary. 
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Appendix A: Monotonicity of the domain-wall 
motion in the retarded ABBM model 

The model defined in (Is]) satisfies Middleton's theorem 
[44| : If the driving is monotonous, wt > Ws for i > s, 
then so is the motion of the domain walF^ 

To prove this, note first that the velocity iit is contin- 
uous. Hence, before any negative velocities can occur, 
there would be an instant t > where iit becomes zero 
for the first time. Then, using ((5|, we get 

rjdtUt — UtF'{ut) + m'^[w{t) — u{t)] — au(t) 

-a I dsf{t - s)u{s) 

= m^w{t)- a I dsf'{t~s)u{s)>0. 

-' — oo 

The first term is w{t) > by monotonicity of the driving 
w. The second term is > since a > 0, f'{t — s) < and 
Us > for all s < t. Similarly one checks that if ilt — 0, 
then dtUt > 0. In other words, at any time where ut = 0, 
the first non-vanishing derivative of Ut is positive. Thus, 
the domain- wall motion is monotonous. 

Note that the vanishing of iiF'iu) for w = is ensured 
if F{u) is smooth. In the case of the ABBM Brownian 
landscape F'{u) is a white noise. One can then consider 



^* Note this is not true for models where domain-wall inertia is in- 
cluded by adding a second-order derivative in the ABBM equa- 
tion of motion [49]. This makes the retarded ABBM model con- 
sidered here quite special. 



a version smoothed at short scale and take the contin- 
uum limit. Alternatively one sees on the formulation 
uF'[u) -i-> y/iiS^{t) in (20) that this is not a problem. 



Appendix B: Position differences in the 
small-dissipation limit 



As discussed in section fVC] instantaneous velocities iit 
are almost surely in the small-dissipation limit 77 = 0. 
However, the distribution of position differences ut — uq 
has a finite limit. The generating function of position 
differences at constant driving velocity v is given by (|26[) : 



,A(«r-uo) = e-Z'o^ '^* '^"' = e"^^^'^) 



(Bl) 



where z{X,T) — f^.Ut- Ut is given by the instanton 



equations (|24|), (|25j) with sources X{t) = \e{t)e{T- t), 
li{t) = 0. In the dissipation-less limit 77 — 0, (251 and 
( 24 1 give an instanton equation similar to (|4l]) 



(1 + a - 2ut)dtUt = ut~u^ + X[ - Xt 



(B2) 



Here we thus work in the limit t^ "C T,Ty{— Sm/v), i.e. 
we set 77 = 0. We express all times in units of t, i.e. our 
variable t is the variable s in section [VA[ In other words, 
we set T = 1, while space units remain such that Sm = 1- 
For i < 0, the sources vanish and the solution is identical 
to that in section [V A| Thus, the integral over u in that 
region has the closed form ( 43 1 , 



utdt = 2uo- + (1 ^ fl) In (1 ^ 'ui-o- ) 



(B3) 



Similarly, the integral over u in the region < t < T is 



Ut At — 



"0+ 



1 + a - 2u , 

z^ -uau 

u — w' — X 



(B4) 



--(a-l)ln(A+(u- 



(a + 4A - l)tanh 



1)")- 



2m- 1 
V1-4A 



V1-4A 



2m 



0+ 



The relationship between uq+ and ut- is given implicitly 
by 



T = 



ds = 



2a tanh 



1 



2u 



«0 + 



2fi-l 



V1-4A 



du 



+ ln(A + (m- 1)u) 



0+ 



(B5) 

The relationship between Uq+ and Uq- , as well as between 
Urp- and Urp+ = is given by the matching conditions 

[(1 + a)uo+ - ul+] - [(1 + a)uo- ~ ul-] = A 

— [(1 + a)uT- — ^T-] ~ ^^- 
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Altogether, this gives a complete (albeit implicit) solu- 
tion for Jj ut in terms of T, A. This allows plotting the 
generating function (Bl), figure 14 It can also be com- 



pared to the result of the standard ABBM model, 



T 



zo(A,T) = - l-Vr^4A 



(B6) 



While the result (B6) holds in the small-dissipation 
(equivalently large-time) limit Tm <C T, a more general 
result was obtained in |TS] (Section II F) for the standard 
ABBM model in the case where t,„ and T are compara- 
ble. 

In figure [Mj j one observes that the slope 
dx\^ z{T,X) — T is independent of the value of 
a. This is also seen from the explicit solution above: 
The instanton equation ( B2 ) , to linear order in A 



(equivalently to linear order in u) simplifies to 
(1 + a)dtut = Ut + X'f- At. 



(B7) 



Its solution is 

1 r°° 



Ut 



d5e-(^-*)/(i+'^'(A;-A,), (B8) 



and its integral is 



dtuf 



ds{X'^-Xs} = XT. 



(B9) 



Thus, z(T, A) = XT + 0{Xf. This means that the av- 
erage displacement, {ut — uq) = vd\\^^„z{T,X) = vT. 
Of course, this is consistent with the fact that the mean 
velocity of the domain wall is fixed by the harmonic well 
rn?{ut — vt). 

On the other hand, from figure [m] (b) one also sees 
that the curvature df\.z{T,X) decreases with increas- 
ing a. Thus, the fiuctuations {ut — uq)^ are decreased 
by retardation effects. This is in agreement with the in- 
tuition of eddy currents slowing the domain wall down 
when it is moving fast, and pushing it forward when it is 
moving slowly. 



Appendix C: Direct perturbation theory in a 

Instead of looking at the slow-relaxation limit t ^ Tm 
or fast-relaxation limit r ^ t„i discussed in sections fVl 



and VI one may attempt to determine the corrections to 
the stationary velocity distribution for a <^ 1 through a 
direct perturbation expansion in a. For this, we need to 
solve the instanton equations (24 1, (25 1 for X{t) = XS{t), 
/i(t) = 0. We take r fixed but a ^ 1 and work, as in 
most of the article, in the units where r,„ = 1, Sm = 1- 
We set 



u{t) = Uo{t) + aui(t) + 0{af 



h{t) = ahi{t) + 0{a)'^ 
(01) 



where uo{t) is the known instanton for the standard 
ABBM model [25], 



hi{t) for < < is then determined by (25 1, 



hi{t) = 



<s~t)~ 



uo{s) 



^ifl,r-i;r-i + l;e-*/-fl-i 



'^^2F,(l,T-';r-' + l;l-\ 



X 

To obtain Ui{t) for t < 0, we need to solve the linear 
equation 

u[{t) ~ ui{t) - uo{t) + 2ui{t)uo{t) + /ii(t) = , (02) 

with Mi(0) = 0. One simple case is r = 1, where 

A 



hi{t) 



1- A 



e*ln(A + e"'(l-A)) 



(03) 



and 



Mt) 



Ae* 



2(A-l)[A(e*-l) + l]2- 
2(A-lfLb(^) 



e*A 



+ (A-l){-A-2(A-l)Li2(^ 

+ {X-l){t-2)t 

+ 3(A-l)ln(A(e*-l) + l) 



2iln 



1-A 



-H2Atln(e~*(l-A))} 
-|-A^e2*ln(A-(A-l)e-*) 

- A(A - l)e*(41n(A - (A - 1)6^*) - l) 

This solution has the asymptotics 

^ ^ e' (2t - 4e* + e^* + 3) 
lim i2i(t)^ ^ ^ ^ln(-A). 



A— >- — oo 



2(e*-l)' 



We see that ui{t) — ^ — oo as A — ^ — oo, for any fixed i < 0. 
Oonsistent with the discussion in sections |V0| and |VD[ 
we recover the result that there is no S{u) contribution, 
Pu^Q — following a kick, and avalanches never end. 

Instead of the stationary velocity distribution, one can 
also consider the stationary distribution of eddy-current 
pressure as in section VA In contrast to (Oil, the con- 



tribution ho{t) of order O(a)'^ to h does not vanish. To 
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(a) z at fixed A = —1.4, as a function of T. Solid blue line: 

standard ABBM model (B6l. Red (dashed), green (dotted), 

yellow (dot-dashed) lines: ABBM model with retardation 

for a = 1,2.3,5.3. 




(b) z at fixed T = 1, as a func tion of A. Dashed blue line: 
standard ABBM model ( B6 1 . Solid blue line: standard 

ABBM model ( [B6| . Red (dashed), green (dotted), yellow 

(dot-dashed) lines: ABBM model with retardation for 

a= 1,2.3,5.3. 



Figure 14. Generating function for position differences z{X,T) as defined in (Bl I 



order 0(0)" the expressions (24 1 and (25) in dimension- 
less units reduce to 

dfUo — uq + Uq + h = 

rdtho -ho = -nS{t). 



Thus 



hoit) = ^e'/^9i-t), 

T 



and the equation for uq becomes 



(C4) 



It is the instanton equation for the standard ABBM 
model (a = 0), but with a time-dependent source Xt — 
^e*/^6l(-t). This is natural, stating that to 0{af, the 
distribution oi iih{t = 0) is the distribution of the observ- 
able /j Xtiit = 7 /°^ dt e^l'^iit in the pure ABBM model. 
The solution of ( |C4| ) is uo(0 == W^' "^^^^^ 

=T't\)"{t) - Ti}'{t) + /ie*/^V(i) 

p*/2 



+C2J.(2e*/(2-)^/;7^ 



Fixing Ci/c2 using the boundary condition {to(i = 0) = 0, 
we obtain 



dtua{t) — In 



V'(O) 
?/;(— 00) 



In oi^i (T,-/iT). 



Thus, the generating function of the stationary distribu- 
tion of h is 



5M/. = [oFi(r,-/.T)]~"[H-0(a)] 



(C5) 



We can now make contact with the result of section fVAl 
Defining /i :— t^^i Vr :— vt and using that 

lim — In qFi{t, —firT ) = ^ 1 + \/l ^ 4/1^ 

r— >C30 7" 



Eq. ( C5 ) in the limit r — >■ 00 reduces to the a — limit of 



(44 1. Its Laplace inverse is given by (45 1 for a = 0. How 



ever, computing the Laplace inverse of ( C5 1 for general 



T, V seems to be doable only numerically. Likewise, we 
did not manage to obtain simple expressions for ui, hi at 
the next-to- leading order, 0{a). 

In the limit v —?' 0^ one defines the density p{h) — 
9t,P(/i) |.u=o+ a-iid one finds 

hp{h) = -LT^^, oFi(r,.r) ' ^^^^ 

where we defined s :== — /i; for the following suppose s > 
0. For half-integer values of t ( C6 I can be expressed in 



terms of elementary functions, for instance for r = 1/2: 



hp{h) = lt;4 



tanh(-\/2s) 



s— >-/l 



2s 



-z2^^s^hv. 



1 



fc=i 



(2fc-l)3 



00 _. 



fc=i 



where Bi is the Jacobi theta function. For small h it 
diverges as p(/i) ~ ^5=''' "^^^' similar to the result (47) 
found in the limit of large r and fixed (not necessarily 
small) a. 
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Appendix D: Perturbation theory in a for the 
backward equation 



In this section, we compute the 6 functions 



.(00) .(10) -(01) 

bative expansion (|164[) for sub-avalanche durations. Ex- 



and h). ,hj. ,h). used in the pertur- 



panding ( 150 ), ( 151 ) in a and ii, and inserting the ansatz 



(164 1, we obtain the following set of equations 



-d,ur 



+ ^r + (.ir)' = -A<5(t-iO (Dl) 



ra,L(i°) + h 



;(io) 



5{t-ti) 



-dt + l + 2u. 



(00) 



-,(10) 



h['°^ = 



(D2) 
(D3) 



-TdM''^+h^^'^ 


-^r = o 


(D4) 


'-d, + l + 2u[°°^' 


u['''>+uf"^ 


- h["'^ = 


(D5) 


-rdM''^ + T^ 


-4''^=0 


(D6) 



Some numerical solutions, and their comparison to this 
perturbative expansion, can be seen in figure |15| 



-(00) 



( Dl ) shows that uj is given by the pure ABBM solu- 



tion (134 1. Similarly, (D2) shows that hi is also given 



(10) • 



by the expression from the pure ABBM model (161), 



^(10) ^ g-(t,-t)/T 



(D7) 



The term ul can still be computed at order a", by 



solving ( D3 1 with the boundary condition uj. 
solution reads 



(10) 



0. The 



^(10) 



dsi exp < / ds2 



re T 
(r2 - 1) [Ae* - (A + 1) 



l + 2^i(0") 



'{tt-Si)/T 



|„A2(r - l)e*(H2) + 2(A + 1)A (r^ - l) e^ 



+t+ti 



- [2A(A+l)T2+2Ar + T + l] 



=*+v+*i 



+ (A+l)2(r + l)eH2t.| 



(D8) 



^(01) 



(01) 



On the other hand, h^ is obtained by solving (D4| with the boundary condition hl^ = 0. The solution reads 



hr^^-- 



r4 e*i--i(A + l 



-A 



Ae^-*. 



(A + l)(r~l) 



(T-i)if / T - 1 1 e*f"*'A 

e^2i^i 1, ;2 --;--- 

T T A -I- 1 



Mr^ / r-1 1 e* *'A 
e - 2^1 1, ;2 



r' \ + l 



(D9) 



Obtaining the higher-order terms for arbitrary r > is 
now complicated, due to the appearance of hypergeomet- 
ric functions. However, the latter simplify significantly 
in the limit t — ;> 1. From now on, we will consider this 
limit only. We then have 



~,(io) _ 



pt-tt 



o* + *i 



^ j-A^e^* + 4(A + l)Ae* 

2 [Ae* - (A + l)e*f ^ 

-he2*'A [-3A - 2(A + 2)i + 2(A + 2)U - 4] 

-2e2*'(i-ii)} (DIO) 



h\ 



(0 



i)_Ae*2^ 
A + 1 



tf - t + In 



(A + l)e*' - Ae* 



(A + l)e*' - Ae*f 



(DID 



(01 



To obtain uj , we need to solve (D5 1 with the boundary 



„-,(oi 



condition u\. ' = 0. The solution reads 



Uf — asi exp s ; 0*2 



dsy 



1 + 2,iif ) 
{uiT - hi":') 



(D12) 



This integral can be evaluated in closed form and gives 
a lengthy expression in terms of logarithms and diloga- 
rithms. However, the expression ( 165 ) for the avalanche 



W 



-,(01) 



duration only depends on J/dtuJ , and on 
These two terms depend only on T := t{ — ti, and are 
simpler: 
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^(01) 



t{ r rtf 

dsi exp < / ds2 



l + 2^(f)]|(uir-^i°^^) 



-r2(A + lf -T (A^ - 2) - (e^ - 1) A(A + 1) + {2T - 3)(A + 1)^ In 1 - 



e-'A 



2(A + l)(-e^A + A + l)' L " "" ■ " " '"" ' '^ -. -v-.-. v- -,v-.-. -^- ^^^ 

+A(3A + 4) In (-e^A + A + l) + 2(A + lfl.h {^^^ + ^(^ + ^^^ ^^ ( aTt) " ^^^ + ^^'^'' ( ATt) 

{A [TA + e^(T{TX + T - 2) + A + 1) + 2T - A - l] 



(D13) 



dtu. 



(01) 



1 



'* 2(A + l)[(e^-l)A-l] 

+ [2e^A((A + 1) ln(A) - (A + 1) ln(A + 1) + 1) + l] In [-e^A + A + 1 



-2e^A(A + l) 



Li2 



1 



A + 1 



Lio 1 



e'X 

A + 1 



(DM) 



To obtain hi , we need to solve (D6| with the boundary condition h\.^ = 0. Note that the formula (165) for the 



avalanche duration only contains h).. . The result for this value again only depends on T :^ t{ — ij: 



hi 



(11) 



T2(A + 1)2 ((e^ - 1) A - 1) + TA (e^ (A^ - 2) - (A + 1)(3A + 4)) + (e^ - l) A(2A + 3)(A + 1) 



(A + l)2(-e^A + A + l) 



(2A + 3) In (-e^A + A + l) 



2Tln 1 - 



e^A \ „. . / e^A 



A + 1 



2Li2 



A + 1 



-2Li2 



A+1 



(D15) 



;(0l 



Finally, /ij. is obtained from (Dll) as 



h'i''^ = 



1 + A 



[T - In (1 + A - Ae^)] 



(D16) 



In figure 15 we show that these perturbative expressions compare well to a direct numerical solution of the backward- 
instanton equations. 



In order to apply (165) we need to compute the leading behaviour as A — > — oo. We get 

T 



lim h. 

A— >- — CO 



(11) 



lim 

A— ^ — CXD 



dt 



^,(01) _ 



.(01) 



-^—^ - Y^ [2Li2 (1 - e^) + (T + l)(r + 2)] 

e^T2 + 2e^Li2 (l - e^) + T + e^ - 1 
2(e^-l) 



2(e^-l)' 



[T^ + 2Li2 (1 - p/) + T 



1] 



lim (-A)"''' exp {-ah,h^^^^\ = (l - g-^)""''' 



(D17) 
(D18) 

(D19) 
(D20) 



Inserting these results into (1651, we obtain (167) 



Appendix E: Avalanche sizes at finite driving 
velocity v in the standard ABBM model 



fixed X ~ vjvm). Let us now consider the stochastic pro- 
cess m(m), i.e. the domain-wall velocity as a function of 
its position. Its first-passage distribution was derived in 
[45] (equivalent to formula (E30) thereP^ and later in 



Similar to the computation of avalanche durations in 
the standard ABBM model at finite driving velocity in 
section fyill A 2[ we can also compute the avalanche sizes. 
The discussion follows closely Ref. [IS] where some ex- 
plicit expressions were obtained in the small-77i limit (at 



1^ V and v' must be replaced by v in the first argument of the 
hypergeometric function there. 
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Figure 15. (Color online) Numerical solution and perturbative 
expansion of ( |150[), ( |151[ ) for a — 0.2, r — 1. Boundary 
conditions are (160 1 with ti = 0, ff = 3.5, /i = 0, A = —1.5. 



Thick (yellow) line: fit from numerical solution. Thin (red) 
line: ht/a from numerical solution. Grey (dashed) line: pure 
ABBM instanton fi*"°\ given by ( |134[ ). Blue (dot-dashed) 

~(01) 



line: h\ given byjDTTJ . Black (dotted) line: u\ + au 



given by (134 1, (D12l 



dMe^"P^(u|uo) = e-5(i-^^^^)(''"-"); 



(El) 



u 


K- 


^ 1 7' ,^1 


- 4Auo 


x/1-4aJ '^'^ ^ 




.5('" 




- 


u 


1 1 7' xT~ 


-4Au 


V1-4aJ'^'^^ 



As w ^- 0, this has a finite limit for u < 1 (the first- 
passage distribution at ?i = 0), 



dMe^"Po(w|wo) 



e-^(i-v^i=4A)^ 



u 


V ( 1 


i ^ „, ./I ^\„-,J 


[2^^ vi-4Ay'^'^^ ^-'-uj 


r(i - v)/T 


1 " fi 1 1 V 


L^ 2 ^^ 1 V1-4aJJ 



Integrating this over the stationary distribution for uq (u) , 
P{uo) = r(t)+i) ^o^~"°' ^^ S6t the Laplace-transformed 
cumulative distribution of avalanche sizes in the form 



.An PJS > U) 

/o°° dS'S'P{S') 



du 



XvB 



(1-4A)^ 

(E2) 



The Laplace transform of the avalanche density p{S) — 
PiS)/[J^ dS"S"P(S")] is obtained by integration by part 
of dEl, 



(1-4A)- 



dS{e^^ -l)piS) = 

vB 



-'!( 



yr^^iA 



(E3) 
A nontrivial check of (|E3|) is that it satisfies the normal- 



ization condition J dS Sp{S) = 1 automatically. 

Using standard relations for the Beta function [S3] , one 
can rewrite (|E3|) as 



' dSie^"^ - l)p{S) ^^(1 _ 4A)-5 fl + ^/T^Ia) 



X B 



V f 1 



(E4) 



Using the substitution r := 1 + Vl — 4A, we can write 
the inverse Laplace transform in the compact form 



p{S) 



ro+ico ^^ sin(7rw) 



27r? 47r 



X re 



\(r-2)rS 



{r-iy-^Biv, 



2-2r 



where tq > 1. 

For the case u = 0, (E3 1 reduces to the known expres- 
sion [m HHj 

/ d5 (e^^ - l)p(5) = - (1 - Vr^4A 

which leads to the standard ABBM avalanche-size den- 
sity p{S) = ^^^e-«/4. 

For any < ii < 1, we can obtain the behaviour of p{S) 
at small S from the A -^ — cx) limit of (E4). In this limit. 



the Beta function tends to a constant and we obtain 

d5 (e^^ - l)p(S') = 

sinTTu, , 1-,. r ,, ,^ 1 /ol / v\ 

= ^^(-4A)^ [1 + 0(|A|)-i/2j B [v,--) 

Inverting the Laplace transform, we see that near 5* = 



p{S) = 5-(3-)/2 



(u - l)r(-z;/2) sin(7r7;) 
47r3/2 



0{S) 



1/2 



(E5) 



Here, B{x,y) is the usual Beta function. (E2) is valid 



since the probability that u = belongs to an avalanche 
of size S is SP[S)/ J^dS'S'P[S') and the conditional 
distribution of u is the n P{ u\S) — 9{S — u)/S. Putting 
this together produces ( |E2[ ). One can check that taking 
the large-A limit and Laplace inverting one recovers ex- 
actly the formula below (E.28) in [JS], valid in the small- 
m limit (at fixed x = v/vm)- 



This is in agreement with previous results [24l ESJ [55] , 
and extends them by giving the prefactor of the power 
law. 



Appendix F: Fixing initial conditions instead of final 
conditions 



In this section, we discuss how to transform (156), a 



formula with a fixed value of hf and a Laplace transform 
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taken with respect to hi, into ( 159 ) , a formula with a fixed to hf 



value of hi and a Laplace transform taken with respect 



We start by integrating (156) over hf, 



dhi e^"'+'''''7'abs(wf , hf, U\ui, hi; U) = 



d/if / dui 



exp<; / [{2 - Ws)us + {I + a + T^'^)] ds + itfiit, \ (271)6 iTh{tf) 



We now invert the Laplace transform from /i = rhiti) to hi using a complex integral, 



oo 

dhi 

-oo Jo 

■ioo 



diii e^"'7'abs(wf , hf, if |ui, hf, tj) 



d/^ exp 



[(2 — Ws)us + (1 + a + r "'^)] ds + ufUt,. > S{hti) exp I —rht-hi j 



exp 



[(2 — Ws)us + (1 + a + T "'")] ds + UfWtj — rhiht. 



^kt{)\hiti: 



=o^(^i) 



(Fl) 



This now is Eq. (159 1. Due to the (5- function, the only 



value of h{ti) that contributes is the one which leads to 
h{tf) — 0. So, the effect of going from a fixed hf to a 
fixed hi in the propagator is a change in the boundary 
conditions for the pair of backward instanton equations 
(1501, (151). When integrating over all hf, we have to 



impose the boundary conditions ( 160 ) 

h{tf) = 0, u(ti) = A. 



In addition, we get the "Jacobian" factor in (159). In the 
pure ABBM case a = 0, as discussed in section |VIIIE1[ 
it just cancels the {tf — ti)/T factor in the exponential. 
For a > it is more complicated. 



Appendix G: Some exact relations for the 
propagator 

Finding general solutions to the forward instanton 

equations (24), ([25|) and the backward instanton equa- 

difficult. However, simple particular 



tions (150) 



( [24|, ([25| ai 
' (|15l| isdif 



solutions exist, where the instanton is constant in time. 
These imply exact relations on particular observables in 
the ABBM model with retardation, which we discuss be- 
low. 



Forward instanton 



A particular solution of (24), (25) is 



Xt^5{t-tf), 
ut = e(tf-t). 



fit = arSit- tf), 
ht^ae{tf-t). (Gl) 



To see the significance of this solution, consider starting 
at fixed Ui,hi at t = 0, and driving with a constant ve- 
locity Wt = V ioT < t < tf. Using (26), and accounting 
for the initial condition as in 



Eq. (4), we get 



g«(tt)+ar/i(tf) _g"i"o+T'ii'io+/o'' dtwtut 



^Ui-\-aThi-\-vt{ 



(G2) 



This implies the following exact relation on the propa- 
gator of the ABBM model with retardation at constant 
driving velocity v, 

/•oo /-oo 

/ diif / d/ife"f+°^'"7'(wf,/if;if|wi,/ii;0) 
Jo Jo 

_ Ui+arhi+vti (C^) 

For the case a = (pure ABBM model) , this relation can 
be checked explicitly using the formula for the ABBM 
propagator, Eq. (19) in [46]. It generalizes similarly to 
the case of arbitrary driving w. 



2. Backward instanton 

Now let us apply the same idea to the backward instan- 
ton, used in section [VIII[ in order to perform calculations 
with an absorbing boundary at u = 0. A particular so- 



lution of (150), (151) is 



Xt^-Sit-tf), 
ut^-e{tf~t), 



Ht = -aTS{t- tf), 
ht^~ae{tf-t). (G4) 



To see the significance of this solution, consider the den- 
sity to arrive at Uf,hf at t = tf, while driving with a 
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constant velocity wt — v for < i < tf, as a function of 



the initial condition at i = 0. Using (154), we get 



M(0)-arh(0) _„«f iitf +-r/if /itj +/(,* dt (2-«;t)fit + (l+a+T ^) 



_ — -iif — ar/if + (i>— 1+a+r )tf 



(G5) 



This implies the following exact relation on the propa- 
gator of the ABBM model with retardation at constant 



driving velocity w, with an absorbing boundary at u = 0: 

diii I d/iie-"'-'^''-7'ate(Mf,/if,tf|Mi,/ii,0) 

(G6) 



— Ui— a/ii + (v — l+a)if 



Again, for the pure ABBM model a = this relation 
can easily be checked using the expression (136) for the 



propagator with an absorbing boundary at u = 0. 
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